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本 发 明 属于 飞机 仿真 技术 领域 ,具体 为 针对 
弹性 飞机 阵风 响应 的 仿真 分 析 方法 和 系统 .本 发 
明 根 据 弹性 飞机 在 阵风 激励 下 的 动力 学 及 结构 
相应 ,首先 根据 飞行 状态 ,对 飞机 六 自由 度 系统 、 
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进行 修正 ,使 之 更 加 符合 弹性 飞机 飞行 的 真实 状 


= 态 ,最 后 ,将 以 上 各 个 模块 连接 成 系统 ,实现 对 弹 
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性 飞机 阵风 激励 下 的 相应 的 分 析 。 该 系统 在 某 通 
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1 .一 种 弹性 飞机 阵风 响应 的 仿真 分 析 方法 ,其 特征 在 于 ,县 体 步 又 包括 : (一 ) 对 飞机 本 


体 飞 行 状态 建 模 , (二 ) 对 飞机 飞行 环境 建 模 , (三 ) 对 飞机 所 受气 动力 以 及 飞机 气动 力 的 弹 


性 影响 建 模 ; (四 ) 将 这 
飞机 的 阵风 相应 ; 


三 个 模型 连 


(一 ) 对 飞机 本 体 飞行 状态 建 模 ,具体 流程 为 : 


飞机 本 体 飞 行 状 态 为 六 自 


机 飞行 的 运动 坐标 系 ; 
HR (S) ` 气流 坐标 轴 
地 面 4 

KRR KELE 


Al £ 


iT 


机 运动 4 


£ 


4 


FEA PR 


AOR 


作 


人 | 


E 标 各 


JE 


并 指向 地 心 ,y 和 


面 内 并 平行 于 飞机 的 设计 各 


IR (S): 
E 标 系 的 原点 为 地 面 上 任 选 的 一 点 (0.) ,其 它 4 
[的 指向 定义 为 x 和 


WARS: F 


ME 


E 标 第 


成 通路 ,建立 模拟 系统 ,最 终 在 Simulink 平 台 上 得 到 弹性 


昌 度 运动 状态 ,为 了 对 飞机 本 体 飞 行 状态 建 模 ,首先 定义 改 
标 系 用 到 以 下 几 种 多 


E 标 


I£ (S) 、 机 体 多 


人 | 


BRE KHL 


E 标 系 按照 欧美 惯用 形式 定义 ,其 中 ， 
标 系 的 原点 者 


质心 (0) 处 且 


在 水 平面 内 指向 某 一 方向 ,z À 


| 垂直 于 地 面 


E 


也 在 水 平面 内 并 垂直 于 x 和 


£ 


站 向 按照 右手 定 则 确 


AE ; x 3 


在 飞机 对 称 平 


线 且 指向 机 头 ,y 


飞机 对 称 平面 且 与 x 3 


x 轴 垂 直 并 指向 机 身 下 方 ,y 3 
对 称 平面 内 的 投影 重合 一 致 ,z Ñ 
体 轴 y, 重合 一 致 ; 
几 个 坐标 系 的 转换 矩阵 如 下 : 
[cosacos sng 0 
体 轴 到 风 轴 :So=| -sinp cosp i 
| 0 0 1 
cos Ocos y 
地 轴 到 体 轴 : S, | 
cos y cos Z 
地 轴 到 风 轴 : Sa ] 


垂直 并 指向 机 身 下 方 :x 和 
垂直 于 0x,z, 平 面 并 指向 机 身 右 方 ;x 
在 飞机 对 称 平面 内 与 x 和 


cos ysin z 


EET KAEKA EN z 3 
上 与 飞行 速度 V 重 合 一 致 ,z 对 称 平面 内 与 


在 


与 飞行 速度 V 在 飞机 


! 垂 直 并 指向 机 身 下 方 ,y $ 


| 与 机 


(1) 


cos Osin y 


(sin øsin Ocosy —cosøsiny) (sm øsin Osin y +cosøcosw) sin gcos 0 
(cos øsin Gcosw +sin øsiny) (cosàsi Osin y -sin øcosy) cosøcos 0 
-sin y 
sin ycos ysin u -sin xcosu sin ysin ysin w+cosY%YcosA cos ysin u 


sin y cos y cos u +sin ysin uy sin ysin y cos 4 — cos cosg COS y cos H 


-sin 0 


o 


(3) 


其 中 ,a 迎 角 ;8B 侧 滑 角 ;96 信 仰角; 了 偏 航 角 ; 由 滚 转角 ; Y 航 迹 倾斜 角 ;x 航 迹 方位 角 ;n 航 


迹 滚 转角 ; 


在 坐标 系 中 对 飞机 质量 特性 建 模 : 


飞机 操 稳 估算 中 的 质量 
生 矩 II、I: 以 及 相对 于 任意 两 个 
性 时 作出 基本 假设 ,飞机 相对 于 对 称 面 几何 外 形 对 称 而 且 内 部 质量 分 布 也 对 称 ,因此 I,、 


轴 的 惯 


I =0; 
yz 


中 ,相对 于 任意 两 个 匀 


转动 惯量 采用 两 种 多 
KRH 


w hR4 


E 标 系 给 出 BUEH 


的 惯性 积 均 为 0， 


E$ 


坐标 系 和 中 心 惯 


Á 


yz 
` 


ER EERI FRS sr 23% , B pe #HXJ-T KHU 
的 惯性 积 1 


XZ 


EER 


] 4| 


RZN 


KER 


\I ;在 研究 飞机 的 运动 特 


Á 


E 标 系 ;在 主轴 坐标 系 


4 标 系 与 体 和 


| 


HiJx 8 Ej z 3: 


21I 


mE x 
2 "E I, 


e e E 


) 


(4) 


4 标 系 的 y 革 
| 方位 角 , 用 a 表示 ,这 个 角度 由 


重合 ,两 坐标 系 
公式 (1) 求 出 : 


CN 110309579 B 


E KOLESI IR E J h RA ABA HAJ KOLE i PE E 3 


权 利 


响 ; 由 一 组 估算 公式 (5) 和 (6) 组 成 : 
ip -WEY 
Lor W(k L,Y 


L, =W(k (b +L,)/27 

Ip =W(b°/78+ H° /33) 
L, =W(L; /29+ H’ /33) 
L, =W(L; /29+b°/78) 
公式 (5) 和 (6) P, T Lp 1, 242 KILE D EER 
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系 的 ,并且 不 考虑 机 落架 的 影 


(6) 


(5) 


系 的 惯性 矩 ,单位 kg-m2;W 为 飞 


MEE, Hkg, bs KILLE RRK B ym, IDE ERE WEK , m, LA KWL 
身 总 长 ,单位 m;k,、k,、k, 为 统计 系数 ,对 于 W 筷 20t 的 飞机 有 k, 二 0.10~0.12,k, 二 0.18~ 
0.19; =0: 26; 3] FW520tBy KAk =0;132e0.45,k =0: 180.19; k, 50:27 AR 


EH FRIZE KALTA ERER, AR O EH TRR EE KALAS E PERE h f ; 
方位 角 的 条 件 下 ,根据 飞机 中 心 惯性 主攻 


EMET KILER 


RE KOUE HR P WEERA PEA; 


_ 3 = 2 
I =l,Sosa+1 sm a 


I, =1, 


L =1 „cos a+] sina 
L. =((ü, -1L,)/2)sin 2a 


建立 六 自 昌 


运动 ,又 有 刚体 的 牵连 


其 中 ,m 为 飞机 质量 ; * 为 速度 矢量 ; FJ Er, 


日 度 运 动 的 动力 学 方程 : 
把 大 地 视 为 平面 ,把 飞机 看 成 刚体 且 有 纵向 对 称 平面 ,飞机 在 运动 中 , 既 有 质点 的 相对 


(8) 


经 推导 得 到 飞机 的 基本 运动 方程 : 


dy 
m A 
(了 


dy 
m —> +ø, v, —@_ v,)= Y 
dt 


dv, 
M aT @,V, -@,v,) = 


+@,v, —@,v,) = X 


(10) 


(11) 


(12) 


MM 为 力矩 矢量 


(7) 


运动 ,根据 牛顿 第 二 定律 ,得 到 力 和 力矩 的 两 个 运动 方程 : 


Ei 


; 了 为 转动 矢 


HE e 
三 三 全 


系 的 惯性 矩 ,利用 公式 (7) 
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do 
六 (13) 
xd “zy yz wd x Z `x 
do do 
Z — + -J 0 o -J C e o J-M 
y dt £ F Ba W ya y (14) 


do 
J —+(J -J )o @ -J_ 0-0)=M 
z dt y X X y xy a z 


RE, vov ov Ax” y zA ARE; 

X、Y、Z 分 别 为 x、y、z 方 向 的 合力 ; 

@ ` o 9, 为 Xx.y、z 方 向 的 角速度 ; 

ML、M、M, 为 x\y、z 方 向 的 力矩 ;飞机 的 转动 惯量 矩阵 为 : 
de dy d 


XZ 


Jald £ Aa 


yx 


J. J J 


在 飞行 载荷 计算 中 ,假设 飞机 的 飞行 速度 V 为 第 数 ,V,sV, 阻 力 等 于 推力 ,得 到 如 下 五 自 
由 度 方 程 : 
dv 


1 = v,@x—v,@, +(Qsina cos P + Ycosa + Zsinasin p) / m — gcos9siny (16) 


Ç: = v,@, —v,@, +(Qsinp — Zcos D) /m+ gcosĝsiny 


do, 5 B 0,0, = B „0,0, + AJ, 2M, + AJ > M, 


do 
m = B 0,0, 


(15) 


(17) 


-B oo +AJ 2M +AJ >M, 


(19) 


&ë J “  % h (20) 


JoJ- J: + J) 
Ih- Jš 


根据 几何 关系 得 到 如 下 三 个 辅助 方程 : 


Bw == 


权 利 
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Li @, Sin y +ø, cos 
+ y SIN y + @, COS y 
< =(@, cosy —@,siny)/ cos9 


一 = @ —tIg9(@j cosy 一 0 sin y) 


将 气动 力 系数 代入 上 述 方程 得 到 如 下 方程 : 
da - o, -B o, — g (WV)Y(Q,z,ó,,M)+ g /Vcos8cosy 


dt 


A + @ cos 
4 „SINY + @, Y 


dt 


r. =@ -tgw cosy — .siny) 


其 中 : 


£ = 


B= 


飞机 姿态 角 与 速度 之 间 的 关系 如 下 : 
v= J: +v +W 


v 
a= E) 


A= sin! (=) 
y 


a 
pei 
0 = sin '(— 
P 


(二 ) 对 飞机 飞行 环境 建 模 


要 求 书 


A Oy +A Q, +g /(WV)[Z,, (9,ó,.e, M)+Z,(Q., B.G, M)+ g /VecosOcosy i 
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(21) 


(22) 


(23) 


(24) 


(25) 


26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) ; 


环境 模型 包括 切 变 风 模 型 、. 素 流 模型 离散 突 风 模 型 ,大 气 参数 模型 和 与 飞行 高 度 相 关 


的 重力 加 速度 模型 ; 
(1) 建立 风 切 变 模型 


风 切 变 的 概念 : 风 切 变 指 的 是 风 矢 量 在 竖 直 方向 上 的 变化 ,以 及 飞机 相对 与 风 矢 量 


其 变化 的 各 种 情况 , 按 航 迹 , 风 切 变 分 为 下 列 4 种 表现 形式 : 
顺风 切 变 ,是 指 水 平 风 的 变量 对 飞机 来 说 是 顺风 ; 
逆风 切 变 ,是 指 水 平 风 的 变量 对 飞机 来 说 的 逆风 ; 
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侧 风 切 变 ,是 指 飞 机 从 一 种 侧 风 或 无 风 状 态 进 入 另 一 种 明显 不 同 侧 风 状态 , 它 有 左右 
之 分 ,使 飞机 发 生 侧 滑 、 深 转 或 偏转 ; 
垂直 切 变 ,是 指 飞 机 从 无 明显 的 升降 气流 区 进入 强烈 的 升降 气流 区 的 情况 ; 


风 切 变 的 天 气 影 响 因素 ,包括 :雷暴 ,锋面 , 超 低 空 急 流 , 低空 逆 温 层 ; 
风 切 变 的 数学 描述 : 
以 地 面 坐 标 系 下 定义 风 切 变 , 即 认为 : 
L 
V = | Ywa 
Lwa (33) 


EEA ORA EE RE SOA IKPA E upa MURAJ TE s RAA v a AE TE ERA E 
下 降 风 为 正 ; 
风速 梯度 表示 为 : 


Quya Olwg Quya 


Oo Oy 6z， uy 


u u 
x Wy Wz 
Ooy Oyy Oy 
Vy, = =| Ver Yey Yie 
j Ox O Oz 
d Yq d i ip i 
OW OWy¿ OWya Wx Wy Wz 


| Óx; Oy; 0z, (34) 


其 中 ,VW 为 纵向 风 、 侧 向 风 和 铅 锤 风 的 水 平 梯度 ,oy ,Vy,, wy, 为 纵 回 风 、 侧 向 风 
和 铅 锤 风 的 侧 向 梯度 ,uy, ,vy, ,Ww 为 纵向 风 、 侧 向 风 和 铅 锤 风 的 重 直 梯度 ; 
建立 风 切 变 模 型 : 

风 切 变 模 型 是 将 风 矢 量 的 变化 加 入 到 系统 仿真 的 模型 中 :基于 Military 
Specification MIL-F-8785C 中 提出 的 算法 ,认为 风速 分 量 的 大 小 仅仅 是 高 度 h 的 函数 , 纵 
向 水 平 梯度 ul ywwww 和 侧 癌 水 平 梯度 uvwww 均 为 零 , 即 : 
_ Auya 
— 

) Vz = Avra 
Az; 

_ Awya 

E 


使 用 如 下 的 方程 来 表示 切 变 风 的 大 小 : 


n) 
u =W Zo ， 3ft< h<1000/ft (36) 


uy; 


Wiz 


(35) 


其 中 ,u 是 切 变 风 在 各 个 方向 下 的 平均 风速 ,W, 是 离 地 6m (20ft) 的 测量 风速 ,h 是 海拔 
高 度 ,z0 是 一 个 常量 ,C 种 飞行 阶段 时 , 值 为 0.0612m; 其 他 飞行 阶段 时 , 值 为 0.816m; 
这 样 得 到 的 切 变 风 是 在 地 面 坐 标 系 下 的 值 ,与 六 自由 度 非 线性 运动 模块 得 到 的 DCM 算 


C 
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阵 相 乘 ,得 到 体 轴 下 的 风速 代入 仿真 模型 


(2) 建立 率 流 模型 ,为 Dryden 连 续 率 流 模型 
大 气 半 流 是 一 个 随机 函数 , 它 与 时 间 和 位 置 有 着 紧密 的 关系 ,这 种 函数 关系 是 基于 大 


量 的 测量 和 统计 数据 来 进行 构建 的 ;假设 : 妆 流 是 平稳 的 、 均 匀 的 ;在 这 种 假设 条 件 下 可 以 


= 


向 


满足 对 飞机 品质 特性 分 析 的 需要 ; 


Dryden 率 流 模型 通过 有 限 差 分 方程 的 滤波 器 ,使 用 有 界 的 白 噪 声 ,把 素 流 对 系统 的 影 


施加 进 仿真 过 程 ， 


大 气 素 流 模型 ,其 对 应 的 指数 型 纵向 相关 函数 为 :f (8) 二 e。““; 基 于 测量 和 统计 数据 求 


== 


得 f (&) 后 ,对 得 到 的 结果 通过 Fourier 变 换 , 求 出 Dryden 模 型 的 纵向 和 横向 频谱 函数 如 下 : 


设计 必须 依据 于 理论 频谱 函数 ,这样 才能 够 保证 所 生成 的 大 气 率 流 数据 频谱 特性 的 正确 


P(Q) = o < PER 
z 1+(L, OQ) 
A 
m 1+(L,0) (37) 
中 (Q) = oM Ls : 2 
z 1+(L,9) 


lm ®(Q) = const 
O 0 


lim (O) < Q? 
ms ， (38) 


其 中 , Q 为 空间 频率 ,单位 为 弧度 / 米 ,L ,L,,L, 为 添 流 尺度 ;0 ,oa JJ Wm aaj 5 ; 


H V w u V w 


使 用 成 型 滤波 器 进行 计算 ,生成 大 气 素 流 数据 ;以 伪 随 机 信号 作为 输入 ,成 型 滤波 器 的 


性 ;成 型 滤波 器 就 是 通过 滤波 环节 将 白 噪 声 转 化 为 期 望 的 有 色 噪 声 ; 


根据 公式 : 
l l-1 l 
o OO=y>laUopo +E YRL Goo.. (o) 
k=l kal p 6 (39) 
(i=1,2...,m) 
其 中 , o 为 时 间 频 率 SH : 
D (@)=G(je)°O (e) (40) 


把 给 定 的 输出 频谱 按照 上 式 分 解 ,得 到 成 型 滤波 器 的 传递 函数 G (s) ;低空 扰动 模型 使 


用 的 是 Dryden 率 流 扰动 模型 ,其 空间 域 频谱 表达 式 ,再 根据 方向 上 飞机 分 速度 V 和 频率 的 关 
系 : 


o = Q V (41) 
得 到 模型 的 时 间 域 频谱 表达 式 为 : 
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i 
V “ z l+(L@1/V) 
l 2L, I+12(L. o / VY 
P = 2 (42) 
V " n h+4Loøoiry} 
p 0) -lo Ls lt 
V " z +4 ory} 
对 上 述 各 式 进 行 分 解 ,得 到 产生 给 定 频 谱 @ ,(o) 所 需要 的 成 型 滤波 器 的 传递 函数 G， 
(s) ;对 于 3 个 分 量 , 求 出 所 需 的 传递 函数 如 下 : 


p, (O) = 


K 
G ç u 
(8) Ts+1 
K, =0, 5 T, = 
zV V 
G (s) = K, (Tas +} 
VE 
K, 一 名 h „Ta = 55 Le = — 
V V V 
Gi ) 区 Kas FD 
S. T 1)“ 
K, = Oy = da = 23L, sda = r 
r V E i (43) 


U EERE uJ A ETKA ARU Ait EL PH a == BJ pk yE 9 sa T AT A AE 28 
函数 的 数学 表达 式 ;为 能 直接 在 程序 中 运算 ,将 滤波 器 所 对 应 的 传递 函数 的 形式 转化 为 可 
以 在 程序 中 以 数值 方法 进行 运算 的 递 推 公式 ,对 于 低空 扰动 的 Dryden 率 流 扰动 模型 , 有 如 
下 递 推 公式 : 

u: Yy =K,(xadt+T,y,)/(T,+df) 

viy =K, l; dt +d? )x > -T dts, a +T} +2T ,df)y,, —T2 y, Ir} +T dt+df) (44) 

wi ya =K [T dry ar), -T dx + QTA +2T dÀ)y, a -Thy ATÈ +T. dt + dt’) 

(3) 离散 突 风 模 型 

对 离散 突 风 来 讲 , 突 风 的 垂直 速度 变化 是 确定 的 ,属于 确定 性 突 风 ;通常 离散 突 风 模型 
有 两 种 :简单 点 的 锐 边 突 风 模型 以 及 复杂 点 的 1-cos 离 散 突 风 模 型 ; 

使 用 在 Military Specification MIL-F-8785C 中 提出 的 1-cos 形 状 的 阵风 模型 ,阵风 
单独 作用 于 每 一 个 轴 或 者 同时 作用 与 三 个 轴 ; 指 定 阵风 强度 ,阵风 长 度 , 以 及 阵风 开始 时 
间 ; 离 散 阵风 模型 以 英尺 / 秒 和 米 / 秒 表示 风速 ;离散 的 阵风 可 以 一 次 或 者 多 次 应 用 ,以 评估 
飞机 多 大 风干 扰 的 响应 ;离散 阵风 的 数学 表达 式 为 : 
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(x <0) 
"a Pa rus cosco) (O<x<qd,) (45) 
Ti F. > y m 


其 中 ,V 为 阵风 强度 ,d 为 阵风 长 度 ,V， ,为 作用 在 体 轴 轴 线 方向 上 的 最 终 风 速 ; 

(4) 大 气 参数 模型 

采用 1976 年 美国 COESA 机 构 发 布 的 美国 标准 大 气 参数 ,这 一 参数 仅 与 飞行 器 所 在 高 度 
有 关 ,包含 了 大 气 绝对 温度 大气 压力 ` 大 气 密度 和 声速 ,在 32000 米 以 下 ; 

E 温度 值 是 由 线性 外 推 得 到 ,压力 值 是 由 对 数 外 推 得 到 ,密度 
和 声速 的 值 则 由 理想 气体 方程 得 到 

(5) 重力 加 速度 模型 

采用 W6S-84 坐 标 系 ,其 原点 是 地 球 的 质心 ,空间 直角 坐标 系 的 Z 轴 指向 BIH(1984.0) 定 
义 的 地 极 (CTP) 方 向 , 即 国际 协议 原点 CI0, 它 由 IAU 和 IUGG 共 同 推荐 ;X 轴 指向 BIH 定 义 的 零 
度 子午 面 和 CTP 赤 道 的 交点 ,Y 轴 和 Z,X 轴 构成 右手 坐标 系 ;WGS-84 李 球 采用 国际 大 地 测量 与 
地 球 物 理 联合 会 第 17 届 大 会 测量 常数 推荐 值 ,采用 的 两 个 常用 基本 几何 参数 ; 

WGS-84 坐 标 系 是 一 个 地 固 坐 标 系 ;有 公式 : 

长 半径 :a 二 6378137 土 2 (m) ; 

地 球 引 力 常数 和 地 球 质量 的 乘积 :GM 二 3986005 X103m*s-2 填 0.6X108ms; 

正常 化 二 阶 带 谐 系 数 :C,, 二 -484.16685X10-6 土 1.3X10-9; 

地 球 重力 场 二 阶 带 球 谐 系 数 :J, 二 108263 X 10-8; 

地 球 自转 角速度 : o = 二 7292115X10-1llrads '+0.150x10-11rads '; 

扁 率 f 二 0.003352810664; 

在 WGS84 中 ,根据 经 度 、 纬 度 和 高 度 ,在 地 球 的 等 位 椭 球 上 确定 重力 加 速度 ; 

(三 ) 对 飞机 所 受气 动力 以 及 飞机 气动 力 的 弹性 影响 建 模 

气动 系数 计算 ,是 根据 飞行 器 的 基本 构 型 ,分 别 计算 气流 坐标 系 下 纵向 基本 气动 系数 、 
横 航 向 基本 气动 系数 、 升 降 舵 产生 的 气动 系数 、 副 经 产 生 的 气动 系数 、 方 向 舵 产 生 的 气动 系 
数 、 襟 翼 产 生 的 气动 系数 、 扰 流 板 产生 的 气动 系数 、 动 导数 产生 的 气动 系数 ,然后 计算 各 部 
分 的 气动 系数 之 和 作为 飞机 整体 的 气动 系数 ; 

计算 出 气流 坐标 系 下 的 气动 系数 后 ,根据 气流 角 将 气动 系数 转化 到 机 体 坐标 系 下 ; 

根据 飞行 器 当前 的 动 压 .重心 位 置 计算 气动 力 和 力矩 沿 机 体 轴 的 分 量 ; 

构建 气动 力 系数 数据 库 : 

先 建立 相关 参数 关于 气动 系数 及 气动 导数 的 多 维 数据 库 ,利用 参数 辨识 校 核 风 洞 实验 
数据 ,相关 结果 与 实验 值 的 对 比如 下 :其 中 ,升力 计算 结果 在 实验 点 满足 5% 的 误差 要 求 ; 阻 
力 系 数 满足 误差 范围 要 求 ;除了 在 小 攻 角 0 附近 由 于 绝对 值 很 小 导致 相对 误差 略 有 超出 ; 

在 仿真 运行 时 ,对 所 有 气动 参数 进行 插值 ,进而 直接 得 到 当前 迎 角 、 操 纵 面 下 机 姻 的 气 
动 参数 ,从 而 输入 下 一 模块 进行 力 和 力 抢 的 设计 ,最 后 与 六 自由 度 非 线性 运动 模块 相 耦 合 ， 
构建 信息 反馈 ; 
分 析 地 面 效 应 对 气动 力 的 影响 : 

由 镜 象 涡 理 论 配 合 经 验 修正 的 综合 方法 出 发 ,考虑 到 飞行 力学 中 的 力矩 配 平 关系 , 形 


9 


CN 110309579 B 权 利 要求 Ë 9/11 页 
成 考虑 地 面 效 应 时 全 机 的 附加 升力 、 阻 力 、 迎 角 和 下 洗 角 的 工程 估算 方法 ,在 给 定 升力 系数 
C 和 离 地 高 度 h 后 ,由 尾 涡 引 起 的 迎 角 减少 量 弧 度 和 附着 涡 引 起 的 升力 系数 变化 量 ,具体 


oC 
ôd = - A (46) 
E 
OC, = — -1 (47) 
š aaa Y | 
其 中 ,各 参量 的 计算 公式 为 : 
F 0.788 
o = exp -2.48 20 | (48) 
2 
p= (2) 2 (49) 


jp/c 
r=—=— n 
8z| (WE) +(1/64)| 
其 中 ,h 为 平均 气动 弦 1/4 处 距 地 高 度 ,b 为 辟 展 , 为 平均 气动 弦 长 
1 2 ' 
NEINN 12 (e Cm) (51) 
£ 


Cu RAAE S ERN EPRA TREI IERAC pN 
FIPE FRR EIGENE ,近似 为 


l 2 8 
N'= + 0.00239 D D | (52) 


(50) 


h h 


构建 气动 弹性 模型 : 

在 气动 弹性 的 计算 中 , 舍 去 对 计算 结果 影响 不 大 的 部 分 ,首先 根据 当前 的 气动 参数 , 计 
算 机 可 上 的 气动 载荷 ,再 将 集中 载荷 根据 经 验 公 式 分 布 在 机 辟 表 面 ,从 而 计算 机 豆 的 二 维 
气动 -结构 耦合 的 结果 ,再 将 机 愤 变 形 量 输 回 到 气动 力 模块 ,得 到 当前 状态 下 的 弹性 影响 的 
气动 力 ; 

KESER AM, OWENE SER; 当 其 结构 的 弹性 特性 可 用 工程 染 理 论 表 达 
时 , 则 是 一 维 静 气动 弹性 问题 ;其 定常 气动 力 用 修正 片 条 理论 ,如 考虑 沿 嗓 展 下 洗 流 变化 ， 
则 用 L.Prandt1 和 J.Weissinger 的 升力 线 理 论 ; 

解 一 维 静 气动 弹性 问题 ,用 积分 方程 矩阵 数值 解法 ;根据 机 囊 几 何 模型 ,积分 方程 为 : 

IEP: 
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f êC 
20)= af” Coned Cit cir Res ant ro) 
j= af c”(y, s) Ga ac + Cuc aean -| c”(y, ndn ,mgdn 
i 88 J 
对 后 掠 圳 ,公式 调整 为 : 
o0)=dj [onide pe 5 ¢an- r'o) 


OO 
(54) 


I Í ôC 7 ] 
了 (y) = af in sa T uac a an Ej z n)+ 2 e (vm) enmedn 
0 0 


a(y)=a,(y)+0(y)+ ó 2 (y) 


, w (55) 
C,(y)=C;(y)+ C; O). L z 


C (yn) 和 Ce (y,n) 94E a 的 过 曲 、 扭 转 刚度 有 关 ; 上 面 各 式 中 ， 
d :远方 来 流速 压 g= 了 OP ;6; 操 纵 面 偏转 角 ,rad;0 (y) :y 处 剖面 顺 航向 弹性 变形 角 ， 


rad;C™ (y ,n) :结构 影响 系数 剖面 作用 单位 力矩 9 方向 引起 剖面 的 弹性 变形 角 ;C“* (y , 

n) :结构 影响 系数 ,nn 剖面 作用 单位 力 z 方 向 引起 剖面 的 弹性 变形 角 8 方 向 ;a,c,e,d,mn,， 
g: 分 别 指 n 站 位 的 几何 迎 角 、 缚 长 ` 气动 中 心 至 刚 轴 的 距离 、 辟 段 重心 至 刚 轴 的 距离 ,单位 展 

长 y 方 向 质量 ,z 向 过 载 、 重 力 加 速度 ; Cz,CisGy xc: 气 动力 系数 ,下 标 L 表 示 升 力 ,M 表 示 力 矩 ， 

AC 表 示 对 气动 中 心 ,上 标 z* 表 示 由 刚性 ,e 表 示 由 弹性 变形 引起 的 ; 

定义 新 的 柔 度 影响 系数 : 


1 
C’(y,7)= c” ,7)+ oC (n) 


1 
CI)=C (y,n)+— C” (v.n) ij 
EDUR JEEE E A EE ERR: 


ejay ce + cc a+) 


()=a|F]| e, + |-n long) 


(57) 


(58) 
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其 中 全 | 为 数值 积分 加 权 和 矩阵 ,如 为 等 间距 间隔 ,用 Simpson 或 梯形 法 则 ;用 升力 线 理论 


时 , 展 向 均 非 等 间隔 ,用 Multhopp 求 积 方程 ,把 区 间 (a,b) 划分 出 Multhopp 站 位 , 则 按 
k 


b 
[fC =Z zÓ) s) sin é, (59) 

a n=1 

如 为 对 称 升 力 分 布 ， o 

í m+l 

ROP n /0 ó, )sn ó, +> 12: ] (60) 
0 


即 可 计算 得 到 |@ EBE. 


2 .根据 权利 要 求 1 所 述 的 弹性 飞机 阵风 响应 的 仿真 分 析 方 法 ,其 特征 在 于 : 

(四 ) 将 三 个 模型 连 成 通路 ,建立 模拟 系统 ,最 终 得 到 弹性 飞机 的 阵风 相应 ,具体 流程 为 
基于 所 建立 的 飞机 本 体 、 环 境 和 气动 力 模 型 ,将 以 上 模型 按 模拟 过 程 中 的 物理 规则 连接 起 
来 ,首先 ,在 每 次 模拟 中 ,根据 六 自由 度 模块 得 到 的 飞机 飞行 状态 量 ,经 环境 模块 计算 得 到 
飞机 当地 所 处 的 风速 温度 、 重 力 加 速度 ;将 这 些 物 理 变 量 与 飞机 飞行 状态 量 结合 ,得 到 气 
流 坐 标 系 中 的 飞机 速 压 、 迎 角 、 侧 滑 角 等 信息 ;把 这 些 信息 放 入 气动 力 数 据 库 进行 插值 ,得 

到 飞机 的 气动 系数 ;在 地 面 效 应 和 气动 弹性 模块 的 作用 下 ,气动 系数 的 值得 到 修正 ;最 后 和 
重力 一 起 ,计算 飞机 飞行 的 合力 以 及 合力 矩 , 传 回 六 自由 度 模块 ,计算 下 一 时 刻 的 飞机 飞行 
状态 ,形成 通路 。 

3. 一 种 基于 权利 要 求 1 或 2 所 述 方 法 的 弹性 飞机 阵风 响应 的 仿真 分 析 系 统 ,由 构建 的 三 
个 模型 :飞机 本 体 飞 行 状态 模型 ,飞机 飞行 环境 模型 ,飞机 所 受气 动力 以 及 飞机 气动 力 的 弹 
性 影响 模型 , 连 成 通路 而 构成 ;在 Simulink 平 台 上 得 到 弹性 飞机 的 阵风 相应 。 
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一 种 针对 弹性 飞机 阵风 响应 的 仿真 分 析 方 法 和 系统 
技术 领域 


[0001] ”本 发 明 属于 飞机 仿真 技术 领域 ,具体 涉及 弹性 飞机 阵风 啊 应 的 仿真 分 析 方 法 和 系 
统 。 


背景 技术 

[0002] ”准确 的 飞机 仿真 模型 对 飞机 控制 系统 的 设计 分 析 有 着 重要 的 作用 。 飞 机 实际 在 空 
中 的 运动 很 复杂 ,不 会 像 教 科 书 上 为 了 说 明 问 题 方便 而 简化 的 运动 模式 那么 纯粹 直接 ,会 
受到 气流 、 速 度 . 温度 、 空 气压 缩 性 等 多 种 因素 的 影响 。 仿 真 分 析 是 无 破坏 性 的 、 可 以 个 性 
化 控制 的 .能够 多 次 重复 的 、 不 受 实际 工作 条 件 限制 ,可 以 实现 的 功能 非常 齐全 的 一 种 实 
验 手段 。 利 用 仿真 实 验 , 可 以 加 深 研 究 人 员 对 实际 系统 的 理解 ,快速 找 出 研制 过 程 中 出 现 
的 问题 的 原因 并 提供 解决 方法 。 
[0003] ”仿真 工具 

[0004] ”本 课题 在 对 飞机 燃油 系统 分 析 计 算 的 基础 上 ,将 采用 MATLAB/SIMULINK 软 件 进 行 
建 模仿 真 ,.MATLAB 是 矩阵 实验 室 (Matrix Laboratory) 的 简称 ,是 一 种 用 于 算法 开发 、 数据 
可 视 化 、 数 据 分 析 及 数值 计算 的 高 级 技术 计算 语言 和 交互 式 环境 .MATLAB 作 为 三 大 数学 软 
件 之 一 , 强 项 就 是 其 强大 的 矩阵 计算 及 仿真 能 力 ,MATLAB 的 应 用 范围 非常 广 , 包 括 信 号 和 
图 像 处 理 、 通 信 、 控 制 系统 设计 、 测 试 和 测量 、 财 务 建 模 和 分 析 , 以 及 计算 生物 学 等 众多 应 
用 领域 。 

[0005] ” MATLAB 软件 中 集成 的 Simulink 工 具 包 极 大 地 方便 了 系统 的 仿真 建 模 过 程 ， 
Simulink 提 供 了 一 种 图 形 化 的 仿真 界面 ,可 以 进行 动态 系统 的 建 模 \ 仿 真 以 及 综合 性 能 分 
析 , 能 够 对 线性 和 非 线性 系统 包括 离散 、 连 续 和 混合 系统 进行 分 析 。SIMULINK 提 供 了 一 个 
交互 式 图 形 坏 境 ，SIMULINK 是 一 个 进行 动态 系统 建 模 、 仿 真 和 综合 分 析 的 集成 软件 包 , 它 
可 以 处 理 的 系统 包括 : 线性 、 非 线性 系统 ;离散 、 连 续 及 混合 系统 ; 单 任 务 、 多 任务 离散 事件 
系统 .Simulink 的 开 放 性 非常 高 ,用 户 可 以 将 系统 的 模型 通过 框图 和 连 线 的 形式 表示 出 
来 ,和 使 用 者 有 很 高 的 交互 性 ,允许 用 户 根据 需要 修改 仿真 模块 的 参数 ,而 且 在 Simul ink 
环境 下 可 以 直接 使 用 MATLAB 软 件 中 的 各 种 分 析 工 具 ,为 仿真 模型 的 性 能 分 析 带 来 方便 。 软 
件 可 以 对 模型 仿真 的 结果 进行 分 析 ， 并 可 以 实时 地 显示 在 可 视 化 窗口 中 .Simulink 内 含 大 
量 的 模块 库 ,使 用 户 能 够 快速 便捷 地 建 立 动态 的 系统 模型 ,只 需要 使 用 鼠标 简单 的 拖 放 和 
添加 模块 间 的 连 线 ,就 能 够 建立 基于 实际 系 统 的 非常 复杂 的 仿真 模型 。 一 般 的 仿真 模型 要 
受到 线性 系统 的 限制 ,而 使 用 Simulink 工 具 包 ， 能 够 建立 更 加 贴近 实际 系统 的 非 线 性 模 
型 ,一 个 自 定义 模块 库 ,并 可 针对 特定 应 用 加 以 扩展 ， 可 应 用 于 控制 系统 设计 、 信 号 处 理 和 
通信 及 图 像 处 理 等 众多 领域 .在 欧美 国家 ,很 多 大 公司 将 产品 投入 实际 使 用 之 前 都 会 进行 
仿真 试验 ,他 们 主要 使 用 的 仿真 软件 就 是 SIMULINK。MATLAB 提供 了 自己 的 编译 器 ,全 面 兼 
容 C++ 及 Fortran 了 两 大 语言 。 因 此 MATLAB 成 为 工程 师 、 科 研 工 作者 最 好 的 语言 .最 好 的 工具 。 
[0006] ”为 了 丰富 Simulink 创 建 仿真 模型 的 功能 ,MathWorks 公 司 陆续 开发 了 许多 有 特殊 
功能 的 模块 库 , 使 得 用 户 可 以 非常 方便 地 建立 仿真 模型 或 者 完成 系统 分 析 。 其 中 主要 用 于 
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航空 航天 仿 真 分 析 的 模块 包 Aerospace Blockset 提 供 了 大 量 的 在 Simulink 环 境 中 使 用 的 
航空 航天 模块 ， 涵 盖 了 航空 航天 飞行 器 设计 中 的 主要 应 用 模型 Aerospace Blockset 和 集成 
了 航空 航天 领域 的 通用 的 标准 模块 ,使 用 户 可 以 方便 地 建立 航空 航天 飞行 器 的 仿真 模型 。 
研究 者 可 以 在 仿真 模型 的 基础 上 进行 全 面 的 系统 开发 和 测试 工作 ,实现 航空 航天 飞行 器 
系统 的 概念 设计 .Aerospace Blockset 通 过 RTW 能 够 自动 生成 实时 代码 ,开展 半 物 理 仿真 研 
F Aerospace Blockset 模 块 库 提供 了 航空 航天 动力 学 和 舱 机 模块 、 飞 行 器 的 推进 系统 、 
空 制 仿 真 系统 ,提供 了 三 自由 度 和 六 自由 度 运 动 模型 库 ,可 以 对 固定 质量 或 者 变质 量 的 飞 
行 器 系统 进行 建 模 与 仿真 ,包含 了 地 球 的 重力 场 人 磁场、 大 气 环 境 和 风 场 的 标准 模型 ,能够 
帮助 用 户 在 设计 中 加 入 航空 飞行 模拟 的 关 键 特性 。 重 力 场 模型 为 1984 年 全 球 测 地 系统 重 
力 模型 WGS1984 和 了 M2000 地 球 磁场 模型 ;大 气 模型 包括 1976COESA 和 标准 大 气 模型 ,以 及 
MIL-HDBK-310 和 MIL-STD-210C 规 范 定义 的 大 气 模型 ; 风 模 型 包括 满足 MIL-F-8785C 和 MIL- 
HDBK-1797 的 离散 阵风 模型 .Dryden 率 流 模型 和 风 切 变 模型 。 用 户 可 以 使 用 各 种 单位 转换 
模块 和 坐标 轴 转 换 模块 以 满足 研究 的 需要 ,同时 还 包 含 动画 模块 ,在 可 视 化 窗口 中 可 以 直 
观 显 示 仿 真 模型 的 动画 。 通 过 在 Simulink 环 境 中 建 模 ， 使 得 系统 结构 非常 清晰 。 与 控制 系 
统 工具 箱 Control System Toolbox 一 起 完成 飞行 控制 系统 的 设计 ,利用 软件 中 的 虚拟 现 
实 工具 箱 Virtual Reality Toolbox 可 以 在 三 维 空间 中 实时 地 显 示 飞 行 器 的 飞行 仿真 姿 
太 


Iùro 


[0007] SIMULINK 仿真 步 又 

[0008] ”在 SIMULINK 提 供 的 图 形 用 户 界面 上 ,只 要 进行 鼠标 的 简单 拖拉 操作 就 可 以 构造 出 
复杂 的 仿真 模型 。 它 外 表 以 方块 图 形 呈 现 , 且 采用 分 层 结构 。 从 建 模 角 度 讲 ,这 既 适 用 于 自 
上 而 下 的 设计 流程 ,又 适用 于 自 下 而 上 的 逆 程 设计 。 从 分 析 研 究 角 度 讲 , 这 种 SIMULINK 模 
型 不 仅 能 让 用 户 知道 具体 环节 的 动态 细节 ,而且 能 让 用 户 清 晰 地 了 解 各 器 件 、 各 子 系统 、 
各 系统 间 的 信息 交换 ,掌握 各 部 分 之 间 的 交互 影响 。 
[0009] ”SIMULIXK 的 每 个 模块 对 于 用 户 来 说 都 是 相当 于 一 个 黑匣子 ,用户 只 需要 知道 模块 
的 输入 和 输出 及 模块 功能 即 可 ,而 不 必 管 模块 是 怎么 实现 的 。 因 此 用 户 使 用 SIMULINK 进 行 
建 模 系统 建 模 时 的 任务 就 是 : 
[0010] ” (a) 分析 计 算 所 需 建 模 的 系统 的 各 个 物理 模型 的 功能 及 原理 ,并 找 出 不 同 元 件 之 
间 的 数 据 关系 ; 

[0011] (b) 如 果 必 要 ,根据 任务 要 求 对 所 需 仿真 系统 进行 适当 简化 , 找 出 物理 系统 必须 建 
模 的 模块 和 它们 彼此 的 内 在 逻辑 联系 ; 
[0012] (c) 在 SIMULINK 中 选择 合适 的 模块 并 把 他 们 按照 自己 的 模型 连接 起 来 。 当 
SIMULINK 模 型 库 中 的 模型 不 符合 自己 的 特定 要 求 时 ,用 户 还 可 以 根据 特定 需求 建立 自己 
定义 的 模块 ,并 把 它 封装 成 一 个 “黑匣子 ”; 

[0013] (d) 进行 调试 和 仿真 , 若 仿真 结果 不 满足 要 求 , 用 户 还 可 以 改变 模块 的 相关 参数 并 
再 次 运行 仿真 ,直至 结果 符合 要 求 。 

[0014] ”气动 弹性 分 析 

[0015] ”实际 情况 中 ,飞行 器 的 结构 并 不 是 绝对 刚性 的 ,飞行 器 外 形 结 构 受 气动 力作 用 会 
发 生 一 定 程度 的 弹性 变形 ,这 种 结构 的 变形 会 反作用 于 气动 力 , 引起 气动 力 的 改变 ,从 而 
导致 结构 外 形 进一步 的 弹性 变形 ,这 样 的 结构 变形 与 气动 力 交互 作用 的 现象 称 为 气动 弹 
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性 现象 ,气动 弹性 效 应 是 飞行 器 研究 过 程 中 不 可 忽视 的 重要 部 分 。 对 弹性 体 飞 行 器 的 研究 
涉及 到 空气 动力 作用 与 弹 性 结构 变形 之 间 相 互 影 响 的 问题 ,对 飞行 器 的 安全 与 性 能 起 着 
关键 的 作用 。 随 着 航空 航天 领域 的 发 展 ,气动 弹性 问题 对 飞行 器 的 安全 和 发 展 的 威胁 逐渐 
突显 ,历史 上 就 曾 发 生 多 起 事故 :1903 年 ,Lanley 的 单 避 机 首次 作 有 动力 的 飞行 试验 时 ,就 
昌 于 气动 弹性 效应 而 发 生 了 机 辟 的 断裂 ， 并 最 终 险 毁 ;1912 年 ,英国 的 Handly Page EBI 
在 作战 时 由 于 尾翼 颤 振 而 坠毁 ;20 世 纪 30 年 代 , 英国 的 " 蛾 号 " 与" 铝 号 "飞机 都 因 发 生 颤 振 
而 失事 :40 年代 ,德国 的 V-2 火 箭 兽 遭 遇 赢 振 而 导致 结构 破坏 .弹性 效应 造成 的 严重 损失 使 
得 人 们 逐步 认识 到 结构 设计 环节 考虑 气动 弹性 效应 的 重要 性 .二 十 世纪 中 叶 以 来 ,研究 者 
们 将 气动 弹性 力学 定义 成 为 一 门 独 立 的 科 学 分 支 并 不 断 发 展 ,其 研究 的 结果 虽然 使 弹性 
问题 造成 的 事故 次 数 有 所 减少 ,但 依然 存在 一 些 损失 严重 的 事故 :20 世 纪 60 年 代 , 北 京 航 
空 学 院 研 制 的 无 人 靶 机 , 因 气 动弹 性 造成 机 体 失 稳 而 坠毁 ;2001 年 ,美国 的 高 超声 速 验 证 
机 X-43 因 控制 舵 面 颤 振 而 首 飞 失 败 , 导 致 整个 Hyper-X 计 划 遭 遇 重 创 ,由 此 可 见 ,对 飞行 器 
弹性 问题 的 研究 已 经 成 为 飞行 器 研究 的 重要 环节 。 现 代 飞 行 器 的 性 能 设计 日 益 追 求 高 速 
度 .高 机 动 性 ,结构 设计 致力 于 减 重 增 效 ,而 飞行 速度 不 断 提高 的 同时 也 带 来 了 复杂 的 气 
动 载荷 ,飞行 器 越 来 越 嘻 现 出 轻 结构 、 高 弹性 、 高 频 响 的 控制 系统 特点 ,这 些 都 会 导致 结 
构 、 空 气动 力 和 控制 系统 等 方面 的 强烈 非 线性 及 其 相互 的 强 耦 合作 用 ,因而 气动 弹性 与 控 
制 系统 相 耦 合 的 问题 在 设计 和 飞行 环节 中 广泛 存在 。 由 于 通常 情况 下 气 动弹 性 现象 会 造 
成 设计 失效 、 结 构 破 坏 的 有 害 影 响 ,工程 界 常 将 其 视 为 不 利 因素 ,这 一 矛盾 突 显 了 气动 弹 
性 力学 与 控制 系统 耦合 研究 的 重要 性 和 必要 性 。 研 究 这 类 问题 ,涉及 到 结构 动力 学 、 非 定 
常 空 气动 力学 和 自动 控制 系统 动力 学 等 多 学 科 领 域 的 交叉 ,具有 较 强 的 综合 性 和 复杂 性 。 
三 者 的 关系 可 以 用 气动 伺服 弹性 力学 三 角形 来 表述 ,如 图 1 所 示 。 

[0016] ”目前 ,气动 弹性 研究 的 主要 方法 有 两 种 : 风 洞 实验 和 数值 模拟 , 风 洞 实验 在 揭示 现 
象 . 探 索 新 课题 方面 准确 可 靠 ,因此 , 风 洞 实验 在 空气 动力 学 的 研究 、 各 种 飞行 器 的 设计 研 
制 方面 ， 以 及 其 他 流体 有 关 的 领域 中 ,都 有 广泛 应 用 ,始终 是 气动 弹性 研究 的 主要 分 支 。 气 
动弹 性 实验 大 体 上 可 以 分 为 两 大 类 :一 类 实验 是 用 来 为 气动 弹性 分 析 取 得 原始 数据 ,如 飞 
行 器 结构 的 刚度 试验 、 结 构 地 面 振动 试验 以 及 地 面 伺服 弹性 试验 等 ; 另 一 类 实验 用 来 获取 
数据 分 析 气 动弹 性 特 性 ,如 赢 振 模型 风 动 试验 、 赢 振 飞 行 试验 以 及 拌 振 试验 等 。 然 而 需要 
指出 的 是 ,气动 弹性 问题 交叉 涉及 众多 学 科 ,试验 研究 的 开展 面临 较 大 的 困难 。 例 如 静 气 
动弹 性 实验 中 ,测量 模型 的 刚 度 会 对 气动 特性 产生 一 定 的 影响 ,通常 的 风 洞 实验 采用 的 模 
型 由 于 选材 问题 都 具有 较 大 的 强度 和 刚度 ,而 真实 的 飞行 器 其 刚度 比 模型 低 得 多 。 因 此 ， 
风 动 实验 中 必需 制造 出 一 种 适用 的 弹性 模型 ,该 模型 能 够 准确 模拟 飞行 器 各 部 件 弯曲 和 
扭转 刚度 ,用 它 进行 风 洞 高 动 压 模 拟 实验 ,用 以 测量 模型 刚度 对 气动 特性 的 影响 首先 , 弹 
性 模型 的 建立 本 身 与 实际 飞行 器 相 比 很 难保 证 弹 性 变形 的 模 态 完全 一 致 ;其 次 , 当 发 生 颤 
振 导 致 模型 结构 破坏 时 ,容易 对 风 洞 设备 甚至 实验 人 员 造 成 伤害 ,具有 相当 大 的 危险 性 。 
风 洞 实验 的 基础 原理 是 相似 原理 ,这 就 要 求 风 洞 流 场 与 真 实 飞行 状态 之 间或 能 满足 所 有 
的 相似 准则 ,或 能 保证 两 个 流 场所 对 应 的 所 有 相似 准则 系数 相 等 .通常 的 风 洞 实验 很 难 完 
全 满足 上 述 两 个 条 件 ,最 常见 的 不 满足 相似 准则 的 情况 是 亚 跨 声速 风 洞 的 雷诺 数 不 符 。 以 
波音 737 为 例 ,该 飞机 在 巡航 高 度 九 千 米 以 上 ,以 每 小 时 927 公 里 的 巡 航速 度 飞行 ,雷诺 数 
达 Re 王 2.4X107, 而 在 3 米 亚 声速 风 洞 实验 中 ,风速 每 秒 100 米 的 情况 下 雷诺 数 仅 约 为 Re 二 
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1.4X106, 相 差 较 大 ,实验 结果 则 会 受到 影响 。 

[0017] ”与 此 同时 , 随 着 计算 机 技术 和 数值 计算 方法 的 发 展 ,数值 模拟 方法 的 研究 范围 以 
及 其 结果 的 精确 度 均 得 到 了 突破 性 的 提高 ,因此 计算 气动 弹性 也 成 为 了 气动 弹性 研究 中 
重要 的 发 展 方 向 ,Schuster 和 Liu 等 给 出 了 较为 通用 的 计算 气动 弹性 的 定义 ,将 采用 数值 
分 析 手 段 . 包 含 线性 和 非 线 性 在 内 各 种 精度 层次 的 气动 弹性 研究 均 列 入 了 计算 气动 弹性 
的 范畴 趾 。 依 据 研究 的 问题 同 结构 、 气 动 方 面 的 非 线性 因素 是 否 相关 ,计算 气动 弹性 又 可 
以 分 为 经 典 线性 气动 弹性 分 析 与 比较 复杂 的 非 线 性 气动 弹性 分 析 。 得 益 于 空气 动力 学 的 
快速 发 展 ,计算 气动 弹性 得 到 了 进 一 步 的 发 展 ,CFD 方 法 采用 了 精确 形式 的 流动 控制 方程 ， 
并 且 数 值 格式 、 离 散 方 法 、 网 格 策略 、 并 行 计算 技术 等 也 逐渐 成 熟 ,目前 ,气动 弹性 计算 方 
法 正 不 断 完善 ,现代 飞行 器 设计 过 程 中 可 以 用 商业 软件 MSC.Nastran 或 自行 编写 的 求解 器 
进行 求解 ,将 气动 弹性 变形 的 影响 引入 到 飞行 器 设计 的 方法 和 过 程 之 中 。 提 高 计算 精度 、 
扩展 分 析 能 力 是 气动 弹性 研究 不 懈 追 求 的 目标 , 因 此 , 面 对 愈加 复杂 的 飞行 器 外 型 .更 加 
优越 的 飞行 性 能 ,研究 包括 非 定 常 气 动力 、 结 构 动力 学 和 控制 方程 等 相互 耦合 的 多 场 耦合 
分 析 方 法 ,是 当前 研究 气动 弹性 问题 的 一 项 重要 任务 。 

[0018] ”近年 来 ,计算 机 技术 和 计算 力学 方法 领域 的 快速 发 展 使 得 采用 计算 流体 力学 进行 
气动 弹性 动力 学 仿真 成 为 可 能 ,研究 气动 弹性 是 研究 固体 形变 在 流 场 作用 下 的 各 种 行为 
以 及 固体 形变 对 流 场 影响 这 二 者 相互 作用 的 一 门 学 科 。1946 年 ,英国 人 Collar 用 一 个 力 的 
三 角形 对 气动 弹 性 问题 作 了 形象 的 分 类 ,如 图 2 所 示 , 三 角形 三 个 顶点 分 别 代 表 空 气动 力 、 
弹性 力 和 惯性 力 ， 从 而 区 分 了 各 学 科 的 研究 范畴 ;三 角形 中 心 的 动 气动 弹性 力学 需要 对 这 
三 种 力 进行 耦合 分 析 ， 是 飞行 力学 、 振 动 和 静 气 弹 等 多 学 科 领 域 的 交叉 , 具有 较 强 的 综合 
EMERE”, 

[0019] FAH AAE a EE e EL OA A PJ i E H PL K RAIA, 
英国 的 Handley Page KEN £ Æ MlZ 9 Je yR I EK S, 这 一 事故 使 得 LanchesteT、 
Bairstow 和 Fages A URII KITAA EHR aj i , BL 22 19224E , MAARRE , i sr >F 4887 
可 以 有 效 地 消除 舵 面 颤 振 ,1934 年 ,Theodorsen 成 功 地 得 到 了 低速 翼 型 在 简 谐 振动 情况 下 
的 非 定常 气动 力 精确 解 ,这 一 理论 被 称 为 Theodorsen 理 论 ,对 气动 弹性 力学 而 言 具 有 划 时 
代 的 意义 加。Dugundji 综 合 叙述 了 飞行 器 的 气动 弹性 的 研究 进展 ,并 且 提 出 了 研究 气动 弹 
性 问题 的 四 个 主 要 手段 :1) 有 效 的 数值 和 理论 分 析 ;2) 风 洞 试验 ;3) 地 面 振动 结构 测试 ;4 
真实 条 件 的 飞行 测 R 

[0020] ”在 经 典 的 气动 弹性 问题 中 ,早期 计算 气动 力 采 用 活塞 . 片 条 和 准 定常 理论 等 比较 
简单 的 方 法 ,文献 中 叙述 了 这 一 阶段 国际 上 的 重大 研究 成 果 。 从 二 十 世纪 七 十 年 代 到 九 十 
年 代 , 核 函数 法 、 小 扰动 位 势 理论 和 偶 极 子 格 网 法 开始 应 用 于 亚 音速 、. 跨 音速 和 低 超 音速 
非 定常 气动 力 的 分 析 中 [9 ,这 些 计 算 气 动力 的 方法 随后 均 有 效应 用 在 了 气动 弹性 的 分 析 
上 。 到 了 二 十 世纪 九 十 年 代 的 中 后 期 , 随 着 计算 机 技术 的 快速 发 展 ,研究 者 们 把 Euler 方 程 
和 N-S 方 程 求解 非 定常 气 动力 的 技术 运用 到 气动 弹性 分 析 方面 ,特别 是 跨 音 速 和 大 攻 角 气 
动弹 性 分 析 方 面 ,取得 了 显著 效果 ,促使 了 气动 弹性 分 析 的 飞速 发 展 " 1 。 

[0021] ”研究 气动 弹性 问题 ,需要 对 非 定 常 空气 动力 学 和 结构 动力 学 (CSD) 等 学 科 进行 交 
又 研究 ， 其 中 包括 了 CFD 与 CSD 之 间 的 耦合 ;CFD 计 算 与 刚体 运动 涉及 动 边界 的 耦合 ;以 及 
CFD 计 算 与 结构 、 运 动 三 者 之 间 的 耦合 等 .其 中 ,求解 结构 动力 学 系统 主要 有 以 下 六 种 数值 
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方法 :1) 振 型 AIMA: 2) 直接 积分 法 ;3) 状态 空间 法 ;4) 时 间 有 限 元 法 ;5) 一 阶 常 微分 方程 
组 初 值 问 题 的 数 值 解法 ;6) 复 模 态 法 。 其 中 直接 积分 法 常用 的 有 线性 加 速度 法 、Newmark 方 
法 、 碍 1son 方 法 等 。 在 研究 直接 模 态 个 加 法 过 程 中 , 邢 誉 峰 等 人 利用 DMSM 筑 略 , 分 析 了 等 截 
面 杆 、 梁 的 碰撞 问题 。 文 献 " 指 出 :这 种 方法 可 以 用 来 获得 结构 弹性 碰撞 问题 的 解析 解 ,而 
且 这 种 方法 不 但 可 以 用 来 分 析 平 动 结构 的 碰撞 问题 ,还 可 以 用 来 分 析 机 构 的 各 种 弹性 锁 
定 问题 ;不 但 可 以 用 来 分 析 结 构 的 点 碰撞 问题 ,还 可 以 有 效 分 析 结 构 的 线 、 面 接触 和 碰撞 
等 问题 .关于 CEFD/ZCSD 的 耦合 方式 ， 流 体 和 结构 系统 进行 独立 求解 ,采用 合适 的 交 做 算法 完 
成 两 个 计算 场 之 间 的 耦合 ,降低 了 各 时 间 步 下 积分 的 复杂 度 , 同 时 也 便于 程序 实行 模块 
化 。 当 前 流 固 耦 合 模型 主要 分 为 紧 耦 合 与 松 耦 合 , 采 用 紧 耦 合 能 够 减少 耦合 响应 过 程 的 时 
间 延 到 ,时间 精 度 得 到 了 提高 同时 增 大 了 时 间 步 长 ， 但 计算 量 也 大 大 增加 ;采用 松 耦 合 在 
每 个 物理 时 间 步 进行 一 次 数据 交换 ,计算 量 与 正常 的 非 定 常 流 场 分 析 相 当 。 此 外 ,Farhat 
和 Zee 等 人 [10] , 叶 正 实 和 蒋 跃 文 等 人 5 还 对 二 阶 时 间 精 度 松 耦 合格 式 进行 了 研究 .气动 
弹性 数值 模拟 过 程 中 , 力 与 位 移 将 通过 插值 技术 在 CED 和 CSD 之 间 进 行 传递 ,因此 存在 两 套 
网 格 :结构 网 格 和 和 气动 网 格 。 非 定常 气动 力 与 结构 位 移 等 参数 需要 在 结构 /气动 网 格 上 进 
行 交 换 , 但 是 结构 /气动 网 格 两 者 的 拓扑 方式 并 不 一 致 。 因 此 必须 要 采用 合适 的 插值 方法 
完成 数据 的 交换 ,插值 方法 需要 满足 以 下 四 点 要 求 :1) 确保 插值 的 光 顺 ;2) 能 够 精确 地 转 
换 两 套 网 格 之 间 的 信息 ;3) 在 网 格 较 稀 朴 的 情况 下 也 能 够 完成 插值 ;4) 能 够 适用 于 复杂 外 
形 。 目 前 比较 常用 的 气动 /结构 网 格 间 插值 方法 有 常 体积 转换 法 (constant volume 
tetrahedron ,CVT) 、 无 限 平板 样 条 法 (infinite plate spline, IPS) 和 薄 平 板 样 条 法 (thin 
plate spline,TPS) 等 。 


发 明 内 容 

[0022] ”本 发 明 的 目的 在 于 提供 一 种 计算 量 小 、 念 真 分 析 性 能 优异 的 弹性 飞机 阵风 响应 的 
仿真 分 析 方法 和 系统 。 

[0023] ”本 发 明 提供 的 弹性 飞机 阵风 响应 的 仿真 分 析 方 法 ,包括 : (一 ) 对 飞机 本 体 飞 行 状 
态 建 模 ，( 二 ) 对 飞机 飞行 环境 , (三 ) 对 飞机 所 受气 动力 以 及 飞机 气动 力 的 弹性 影响 建 模 ; 
(四 ) 将 这 三 个 模型 连 成 通路 ,建立 模拟 系统 ,最 终 在 Simulink 平 台 上 得 到 弹性 飞机 的 阵风 
相应 。 

[0024] (一 ) 对 飞机 本 体 飞行 状态 建 模 

[0025] ”飞机 的 气动 数据 和 非 线性 运动 方程 是 飞机 非 线性 动力 学 模型 的 基础 .飞机 的 气动 
力 和 力矩 系数 、 静 稳定 性 和 动 稳定 性 导数 是 由 飞机 的 基本 构 型 和 飞行 状态 确定 的 ,操纵 导 
数 或 操纵 引起 的 气动 系数 变化 与 操纵 面 的 设计 和 操纵 输入 量 密切 相关 。 在 本 节 中 ,由 于 飞 
机 本 体 飞 行 状态 为 六 自由 度 运 动 状态 ,为 了 对 飞机 本 体 飞行 状态 建 模 ,因此 首先 要 定义 飞 
机 飞行 的 运动 坐标 系 ， 并 在 机 体 坐 标 系 中 ,对 飞机 质量 特性 分 布 和 飞机 转动 惯量 进行 建 
模 ,最 后 推导 飞机 六 自由 度 动 力学 方程 ,最 后 在 Simulink 上 得 以 实现 。 

[0026] ”飞机 运动 坐标 系 的 定义 

[0027] ”作用 在 飞机 上 的 各 种 力 和 力矩 的 产生 原因 各 不 相同 ,因此 ,选择 合适 的 坐标 系 来 
方便 确切 的 描述 飞机 的 空间 运动 状态 是 非常 必要 的 。 本 文中 主要 用 到 以 下 几 种 坐标 系 :地 
面 坐标 轴 系 〈$.)、 机 体 坐 标 轴 系 (Su 、 气 流 坐 标 轴 系 (5,) 、 稳 定性 坐标 轴 系 (Ss) 。 坐 标 系 按 
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照 欧美 惯用 形式 定义 ,如 图 3 所 示 。 其 中 ,地 面 坐标 系 的 原点 为 地 面 上 任 选 的 一 点 (00 ,其 
它 坐 标 系 的 原点 都 取 在 飞机 质心 (0) 处 且 坐标 系 与 飞机 固 连 .各 坐标 轴 的 指向 定义 为 xs 划 
在 水 平面 内 指向 某 一 方向 ,z 轴 垂 直 于 地 面 并 指向 地 心 y 轴 也 在 水 平面 内 并 垂直 于 x 3 
(指向 按照 右 手 定 则 确定 ) ;x, 轴 在 飞机 对 称 平面 内 并 平行 于 飞机 的 设计 轴线 且 指 向 机 头 ， 
y Ba LT KIL 对 称 平面 指向 机 身 右 方 ,z 轴 在 飞机 对 称 平面 且 与 x 轴 垂直 并 指向 机 身 下 
方 ;x, 轴 与 飞行 速 度 V 重 合 一 致 ,z 对 称 平面 内 与 x 轴 垂 直 并 指向 机 身 下 方 ,y 轴 垂 直 于 
Ox z 平面 并 指向 机 身 EN x 轴 与 飞行 速度 V 在 飞机 对 称 平面 内 的 投影 重合 一 致 ,z_ 轴 在 
飞机 对 称 平面 内 与 x. 轴 垂 直 并 指向 机 身 下 方 ,y 轴 与 机 体 轴 y 重合 一 致 。 

[0028] “” 几 个 主要 坐标 系 的 转换 矩阵 如 下 : 
lcosgcosB sin £ I 


(1) 


[0029] ” 体 轴 到 风车 


:Su=| -snp cosp 0 
0 0 1 


cosOcosy cos @si y -sin 0 
[0030] “地 轴 到 体 轴 : S, = (sn Asin Gcosy —cosósn iz) (sin ġsin Osin y +cosøcosw) sin ócos0 | (2) 
(cosøsin @cosu +sin Gsinw) (cosàsm sin y -sin cosy) cosøcos 


(3) 


[0031] “地 轴 到 风 旬 


snycosxsin -sin xcosu sin ysin sin 4 +Cos y cos cosysin u 


cosy cos Z cosy sin z 一 Sn y 
] : Sa 5 


sin ycos ycos u +sin ysin uy sin ysin y COS- Cosy Cosy cosycosu 
[0032] FE, aii f: BME fA; OMNA Dm: b 23 8; y 航 迹 倾斜 角 ;x 航 迹 方位 
角 ;u 航 迹 滚 转角 。 

[0033] ”在 坐标 系 中 对 飞机 质量 特性 建 模 

L0034] ”飞机 操 稳 估算 中 的 质量 特性 主要 指 转动 惯量 参数 ,也 就 是 相对 于 飞机 体 轴 系 三 个 
坐标 轴 的 惯性 矩 I, II 以 及 相对 于 任意 两 个 坐标 轴 的 惯性 积 1L. II 如 前 所 述 ,在 研 
究 飞 机 的 运动 特性 时 作出 基本 假设 ,飞机 相对 于 对 称 面 几 何 外 形 对 称 而 且 内 部 质量 分 布 
也 对 称 , 因 此 1 = 1 0. 

[0035] “飞机 设计 中 的 转动 惯量 通常 采用 两 种 坐标 系 给 出 , 即 体 轴 坐 标 系 和 中 心 惯性 主攻 
坐标 系 。 在 主轴 坐标 系 中 ,相对 于 任意 两 个 坐标 轴 的 惯性 积 均 为 0 ,主轴 坐标 系 与 体 轴 坐标 
系 的 y 轴 重 合 , 两 坐标 系 的 x 轴 与 z 轴 相差 一 个 角度 ,这 个 角度 称 为 主轴 方位 角 , 用 a 表示 ,这 
个 角度 可 以 由 公式 (1) 求 出 : 

27, 
J. = 1, 
[0037] ”对 飞机 转动 惯量 的 估算 公式 都 是 相对 飞机 中 心 惯性 主轴 系 的 ,并 且 没 有 考虑 机 落 
架 的 影 响 。 常 用 的 一 组 估算 公式 由 (5) 和 (6) 组 成 。 

L, =W(k bY 
[0038] J, =W(k, L,Y (5) 
L, =W(k.(b+ L.)/2)° 


1 
[0036] a= arctg( ) ¿42 
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L, =W(P'/T8+ H° 433) 
[0039] Z, =W(L; /29+H?° 1/33) (6) 


L, =W(L; /29+b°/78) 

[0040] AR (5) (6) ,1 Lp 1, HAP KOLP CEE AR TERE (单位 kg-m2) ; 
W 为 飞机 重量 (单位 kg) bA KILLE RK (单位 m ,LL 为 不 含 空 速 管 的 飞机 总 长 (单位 m ,L; 
为 飞机 机 身 总 长 Eim) ;k,,k,,k, 为 统计 系数 ,对 于 W 三 20t 的 飞机 有 k, 二 0.10~0.12，k, 二 
0.18~0.19,k, 二 0.26; 对 于 W220t 的 飞机 有 k, 二 0.13~0.15,k, 二 0.18~0.19,k, 二 0.27。 公 
式 〈2) 适用 于 战斗 机 等 小 展 弦 比 飞机 的 惯性 矩 估算 ,公式 (3) 适用 于 中 型 运输 机 、 旅 客机 等 
大 展 弦 比 飞 机 的 惯性 矩 估 算 。 

[0041] 在 确定 了 飞机 主轴 方位 角 的 条 件 下 ,根据 飞机 中 心 惯性 主轴 系 的 惯性 矩 可 以 利用 
ARMOR 出 飞机 体 轴 系 中 的 惯性 矩 和 惯性 积 。 


Š > 2 
I =1„cos a+1„sn“a 


1, = 1 
[0042] Ü” ” f (7) 
1 =l, cos? a+ly sin” a 


L, =((1, -1,)/2)sin 2a 
[0043] ”在 早期 设计 阶段 ,飞机 的 主轴 方位 角 往往 难以 确定 ,对 于 轻型 飞机 或 展 弦 比较 小 
的 战斗 机 通常 可 以 忽略, 即将 主轴 方位 角 假 设 为 0, 此 时 飞机 相对 于 中 心 惯 性 主轴 系 与 机 
体 轴 系 的 转动 惯量 相同 , 量 惯 性 积 的 值 为 0; 对 于 大 型 飞机 ,了 ,的 作用 较 明显 ,必须 考虑 主 
轴 方 位 角 , 初 步 设计 中 通常 取 某 一 近似 值 。 
[0044] ”建立 六 自由 度 运 动 的 动力 学 方程 
[0045] ”把 大 地 视 为 平面 ,把 飞机 看 成 刚体 且 有 纵向 对 称 平面 ,飞机 在 运动 中 , 既 有 质点 的 
相对 运 动 , 又 有 刚体 的 牵连 运动 ,根据 牛顿 第 二 定律 ,得 到 力 和 力矩 的 两 个 运动 方程 : 


[0046] m= =F 

dt (8) 
[0047] MT, 

dt (9) 


[0048] ”其 中 ,m 为 飞机 质量 ;5 为 速度 矢量 ; FAIRE; 下 为 力矩 矢量 ; 工 为 转动 矢量 。 
[0049] ”经 推导 得 到 飞机 的 基本 运动 方程 : 


[0050] a 
dt 


+@,v, —@,v,)= X 


(10) 
dv 
[0051]  m(—+@,.v, —@,v,)=Y 
dt (11) 
dv 
二 十 = =Z 
[0052]  m( a @,V,) D 
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do, do, 
[0053] J muL A -J Ma a -7 €C 7 aa )=M 
X d Z X` dt xZ X (13) 
a dw 
[0054] J 一 -J DO, m. 一 (—*-øo w )= M 
Y dt Jy dt y z y (14) 
do, s P 
[0055] J —<+(J -J -J -0o )=M 
zy oy ry) M, (15) 


[0056] ”其 中 ,vv,、v, 为 x\y、z 方 向 的 速度 ; 

[0057] XY ZAJ xy z AE; 

[0058] o ono xy zA B 8 EE; 

[0059] M.M, M ,为 x\y、z 方 向 的 力矩 ;飞机 的 转动 惯量 矩阵 为 : 


Hl 


Z k a| 
[0060] J=|J, Pi Íy le 
Jx Jy J, 


[0061] 在 飞行 载荷 计算 中 ,假设 飞机 的 飞行 速度 V 为 常数 ,V,sV, 阻 力 等 于 推力 ,得 到 如 
下 五 自由 度 方程 : 


dv, 
[0062] Pi @x— v,@, +(Qsina cos P + Ycosa + Zsinasinp)/ m — gcosQsiny 
(16) 
dv . : 
[0063] =v,@, —v,@, +(Qsinp — ZcosD))/ m+ gcos9siny hi 
do 
[0064] ”一 = 有 OO -B 0,0, +AJ XM, +4J, 2M, 
dt (18) 
do 
[0065] “一 一 = B.,@, o, - B oo, + AJ. 29M, +AJ >M, 
dt (19) 
g J. J 
[0066] 2 F “0 GQ, + x (œ 0)+ EM, 
Z z z (20) 
[0067] 其中. 
l 
A= 
[0068] I5- 
JT? — JJ, + J? 
[0069] B,=—— 
JJ; — J, 
V — Jo gO 
[0070] B,= sidi 
JJ, Jo 
So — JJ: + A 
[0071] Bw = JIo — J: + Jy) 


.无 五 一 .72 
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[0072] 


[0073] 


[0074] 


[0075] 
[0076] 


[0077] 


[0078] 


[0079] 


[0080] 


[0081] 


[0082] 


[0083] 


[0084] 


[0085] 


[0086] 


[0087] 


[0088] 


[0089] 


说 HH + 
根据 几何 关系 得 到 如 下 三 个 辅助 方程 : 


9/20 页 


=@, siny +, cos y 


dt (21) 
ago cosy —@ sin y)/ cos9 

dt à i (22) 
A —tgeQ9(@ , cosy —@, sin y) 

4 š i (233 


将 气动 力 系数 代入 上 述 方 程 得 到 如 下 方程 : 


da - ©, -P o, — g (WV)Y(Q,a,ó,, M) + g/ VeosOcosy b 
“asura cosy 
4 | (25) 
dp 
r i +a o, + g/(WV) Zy (Q, ô, a, M)+Z(0, p,a, M)+ g! VcosOcosy “pas 
a 一 如 OO cosy — ,siny) 
f I i (a7) 
其 中 : 
ale Wa Hra 
Pe, A da F (28) 
V, d dy, 
B= £ p Z V, 
V, dt dV; 
飞机 姿态 角 与 速度 之 间 的 关系 如 下 : 
_ J 25 2 5 
Wal Ve Fy +A (29) 
a= tg C) 
Vy (30) 
“a 
B= sin '(—) 
v (31) 
0 = sin (22 
v (32) 


基于 以 上 方程 ,在 Simulink 中 建立 飞机 的 非 线性 动力 学 模型 ,然后 加 入 输入 输出 


接口 就 可 以 作为 飞机 的 本 体 非 线性 模型 在 飞行 仿真 中 使 用 .这 样 建立 的 模型 还 没有 加 入 


反馈 控 


判 环节, 只 有 从 输入 到 输出 的 单 向 传递 关系 。 而 这 里 非 线性 动力 学 模型 全 部 采用 


Simulink 中 已 经 建立 运算 关系 的 封装 模块 ,联合 其 他 动力 学 及 模块 进行 有 反馈 的 非 线性 


动力 学 仿真 ， 


[0090] 
[0091] 


质量 模块 的 计算 结果 也 作为 输入 参数 参与 其 中 。 
(二 ) 对 飞机 飞行 环境 建 模 
在 建立 飞机 本 体 非 线 性 动力 学 模型 时 ,为 了 提高 仿真 分 析 的 逼真 度 ,在 全 机 仿真 
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模型 中 需 要 建立 较为 真实 的 环境 模型 ,环境 模型 包括 切 变 风 模 型 、. 素 流 模 型 、 离 散 突 风 模 
型 大气 参数 模型 和 与 飞行 高 度 相 关 的 重力 加 速度 模型 。 

[0092] ，” 风 切 变 模块 和 率 流 模块 参见 图 4 所 示 。 

[0093] (D) 风 切 变 模型 

[0094] ” 风 切 变 指 的 是 风 矢 量 在 竖 直 方向 上 的 变化 , 根 飞 机 相对 与 风 矢量 及 其 变化 的 各 种 
情况 , 按 航 迹 可 以 把 风 切 变 分 为 下 列 4 种 表现 形式 : 

[0095] ”顺风 切 变 (TAIL WIND SHEAR) 指 的 是 水 平 风 的 变量 对 飞机 来 说 是 顺风 。 顺 风 切 变 
使 飞 机 空 速 减 小 ,升力 降低 ,飞机 下 沉 ,这 是 一 中 比较 危险 的 风 切 变形 式 。 

[0096] ”逆风 切 变 (HEAD WIND SHAEAR) 指 的 是 水 平 风 的 变量 对 飞机 来 说 的 逆风 。 首 风 切 变 
使 飞 机 空 速 突然 增 大 ,升力 突然 增 大 。 飞 机 抬升 ,危害 相对 轻 些 。 

[0097] ” 侧 风 切 变 (CROSS WIND SHAEAR) 指 的 是 飞机 从 一 种 侧 风 或 无 风 状 态 进 入 另 一 种 明 
显 不 同 侧 风 状 态 , 它 有 左右 之 分 ,使 飞机 发 生 侧 滑 。 滚 转 或 偏转 。 

[0098] “垂直 切 变 (VERTICAL WIND SHEAR) 指 的 是 飞机 从 无 明显 的 升降 气流 区 进入 强烈 的 
升降 气流 区 的 情况 ,特别 的 强烈 的 下 降 气流 ,往往 有 很 强 的 狮 发 性 ,强度 大 ,使 飞机 突然 下 
沉 , 危 险 很 大 ,参见 图 5。 
[0099] 天气 因素 的 影响 : 

[0100] ”雷暴 

[0101] 雷暴 是 产生 强烈 低空 风 切 变 的 重要 天 气 。 雷 暴 单 体 下 方 的 下 蝶 气 流 在 相当 宽 的 范 
围 内 可 以 造成 由 下 击 暴 流 和 雷暴 外 流 组 成 的 两 种 不 同 的 风 切 变 。 一 种 是 雷暴 单 体 中 心 附 
近 下 方 的 下 击 暴 流 切 变 , 他 表现 为 范围 小 ,生命 期 短 ,强度 大 的 特征 。 另 一 种 是 下 击 暴 流 接 
近 地 面 时 转化 为 强 烈 的 冷 空气 外 流 , 它 向 四 外 传播 , 沿 雷暴 单 体 前 进 方向 伸 出 15 一 25km， 
并 使 暖 湿 空气 入 流 抬 升 形成 阵风 锋 ,在 雷暴 下 方 大 范 围 内 引起 180 度 的 风向 变化 ,表现 为 
强 顺 风 切 变 和 强 逆风 切 变 。 由 于 其 中 一 部 分 强风 切 变 区 向 前 伸展 远离 雷暴 主体 ,不 易 察 
觉 , 对 飞行 安全 威胁 很 大 。 当 飞机 穿越 这 类 风 切 变 时 ,先是 逆风 增 大 ,而 后 逆风 对 飞机 速度 
减 小 至 零 ,并 有 强 的 下 降 气流 , 紧 跟 着 又 是 大 顺风 ,这 种 风 切 变 对 飞行 危害 是 最 大 的 ,飞机 
在 这 一 过 程 中 先是 空 速 增 大 ,使 飞机 有 利于 获得 高 ,然后 空 速 迅 速 减 小 ,使 飞机 迎 角 减 小 
也 造成 升力 降低 , 紧 跟 的 大 顺风 会 使 风速 再 次 减 小 ,如 果 空 速 减 到 一 定 程度 会 使 飞机 掉 高 
度 甚至 险 地 。 

[0102] ”锋面 

[0103] ”锋面 是 产生 风 切 变 最 多 的 天 气 系统 ,锋面 两 侧 气 象 要 素 有 很 大 差异 ,锋面 过 渡 区 
的 垂直 结 构 , 是 产生 风 切 变 的 重要 条 件 。 一 般 当 锋面 两 侧 的 温度 差 大 于 等 于 5 度 ,锋面 移动 
速度 大 于 等 于 15m/s 时 ,都 会 在 锋面 附近 产生 对 飞行 有 影响 的 低空 风 切 变 ,其 中 尤其 以 冷 
锋 型 低空 风 切 变 危害 较 明 显 , 但 这 种 低空 风 切 变 一 般 持 续 时 间 较 短 . 在 冷风 后 的 偏 北 大 风 
区 内 往往 也 存在 较 严 重 的 低空 风 切 变 。 如 1973 年 在 波士顿 罗 津 国际 机 场 一 架 坠 毁 的 DC-10 
就 与 锋面 有 关 , 男 外 冷 锋 附近 会 经 常 形成 雷暴 ,也 会 产生 强 的 低空 风 切 变 , 因 此 我 们 不 能 
对 它 忽视 。 

[0104] ” 超 低 空 急流 

[0105] ”在 1500 米 以 下 ,中 心 位 置 常见 于 120--160 米 之 间 的 超 低 空 急流 常 伴随 产生 低空 风 
切 变 ， 最 大 风 切 变 位 于 急流 轴 之 下 ,我 国 曾 观 测 过 到 的 最 低 四 十 七 米 的 超 低 空 风 切 变 。 超 
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低空 急流 出 现时 ,一 定 伴 有 逆 温 , 正 是 逆 温 层 阻挡 了 在 其 上 的 大 尺度 运动 与 地 面 之 间 的 动 


[0106] 
[0107] 


量 交换 ,为 逆 温 层 以 上 动 
层 上 和 急流 的 形成 和 维 FF oH 


上 时 
低空 逆 温 层 


KA 


青 夜 ,在 低空 易 存在 革 


囊 射 逆 温 


层 大 风 癌 下 层 的 传输 ,使 地 面 风 很 弱 , E 
引起 的 低空 
切 变 稍 弱 些 .这 种 风 切 变 存 在 时 往往 逆 ; 
温 层 时 会 造成 上 升 或 损失 高 度 , 如 果 逆 ; 
易 造 成 飞机 险 地 ,或 冲 出 跑道 。 这 种 风 切 变 存 在 时 往往 逆 ; 
大 做 ,飞机 穿越 逆 》 


层 , 并 


量 的 累积 和 储存 提供 了 


HBH ,地 面 月 


了 条 件 ,同时 减缓 了 动量 的 耗 散 ,利于 逆 温 


风 。 


常 伴 有 低空 急流 ， 


由 于 稳定 的 逆 》 


昌 层 阻碍 了 上 


昌 层 上 下 形成 风 的 垂直 切 变 ,单纯 


RSI- 


mja ,下层 风 平 浪 前 


昌 度 时 会 造成 上 升 或 损失 高 度 , 如 果 逆 温 


风 切 变 的 强度 较 雷 暴 和 锋面 引起 的 风 切 变 弱 的 多 ,也 比 超 低 空 急流 引起 的 风 
,而 上 层 狂风 大 做 ,飞机 穿越 逆 
昌 层 很 低 , 在 飞机 着 陆 时 ,飞行员 处 置 不 及 时 ,很 容 


昌 层 , 下 层 风 乎 浪 间 
度 很 低 ,在 飞机 着 陆 时 ,飞行 员 


,而 上 层 狂风 


处 置 不 及 时 ,很 容易 造成 飞 机 附 地 ,或 冲 出 跑道 ,此 外 当 机 场 周围 的 环境 和 
时 ,也 会 产生 对 飞机 的 起 飞 , 着 陆 有 影响 的 低空 


山地 高 差 大 ,水 域 面 积 大 ,机 场 附 近 高 大 建筑 物 群 等 ， 
条 件 下 。 它 出 现 的 高 度 一 般 较 低 , 并 且 通 常 


伴随 高 度 


也 形 e 


风 切 变 。 如 地 形 波 , 较 大 水 陆 界 面 等 
AAPEA JEBE 
兽 加 而 风速 减 小 ,飞机 在 降 过 程 中 容 


易 遭 遇 这 种 风 切 变 , 因 在 起 降 过 程 中 飞机 高 度 低 , 速 度 小 ， 而 我 们 通常 是 季风 起 落 ， 逆风 的 


减 小 使 飞机 空 速 减 小 ,如 果 无 法 使 飞机 速度 加 速 来 应 付 风 对 飞机 空 


测 , 重 则 导致 机 坠 地 ,造成 严重 后 果 。 


[0108] 
[0109] 


[0110] 


[0111] 


风 切 变 的 数学 描述 ; 


为 了 研究 方便 ,这 里 在 地 面 4 


Uya 
Vsa =| Vwa 


Wwa c33) 


每 个 风速 分 量 定 义 为 :水 平 风速 


量 w,， 下 降 风 为 正 。 
[0112] ”风速 梯度 可 以 表示 为 : 
| Buwa OUuwa Quya 
Ox O, 0z， Uy, 
Vwa Vwa Owa 
[0113] Vy, = =| Vx 
4 Oxy Oy 0， 
Owwa OWwg OWwa Wire 
Oy O G, 


[0114] 
AANER 
[0115] 
[0116] 


Specification MIL-F-8785C 中 提出 的 算法 ,这 里 认为 风速 分 量 的 大 小 仅仅 


其 中 ,Vy ;Wi 为 纵向 风 、 侧 向 风 和 铅 包 
锤 风 的 侧 向 梯度 ,uv , Vy Wy 为 纵向 风 、 俱 


建立 风 切 变 模型 : 


E 标 系 下 定义 风 切 变 , 即 认为 : 


uo MAAE: WANE 


(34) 


H 


向 左 为 


E Vyd , 


速 的 影响 , 轻 则 影响 目 


正 ; 铅 锤 风 分 


E 风 的 水 平 梯度 ,Vy ,Wy 为 纵向 风 、 侧 向 
| 向 风 和 铅 锤 风 的 垂直 梯度 。 


风 切 变 模块 就 是 将 风 矢 量 的 变化 加 入 到 系统 仿真 的 模型 中 。 基 于 Military 


数 ,纵向 水 平 梯度 uvwwww 和 侧 向 水 平 梯度 ul ,yw ww 均 为 零 , 即 
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[0117] 


[0118] 


[0119] 


[0120] 


[0121] 
[0122] 


的 DCM 和 矩 FBF4H3F849 301073 


[0123] 
[0124] 


数 , 它 与 时 间 和 位 置 有 着 
构建 的 。 在 航空 工 程 应 用 


(35) 


使 用 如 下 的 方程 来 表示 切 变 风 的 大 小 : 


h 
Zo 3fi<h<1000fi 


明 + 


(36) 
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其 中 ,u, 是 切 变 风 在 各 个 方向 下 的 平均 风速 ,Wo 是 离 地 6m (20ft) 的 测量 风速 ,h 是 
海拔 高 度 ,z0 是 一 个 常量 ,C 种 飞行 阶段 时 , 值 为 0.0612m; 其 他 飞行 阶段 时 , 值 为 0.816m。 

切 变 风 在 竖 直 方向 上 的 分 布 ,参见 图 6 所 示 。 
这 样 得 到 的 切 变 风 是 在 地 面 坐标 系 下 的 值 , 需 要 与 六 自由 度 非 线性 运动 模块 得 到 


(2)Dryden 连 续 率 流 模型 


I 下 的 风速 代入 仿真 模型 。 


实际 的 大 气 率 流 是 十 分 复杂 的 物理 现象 .为 了 使 飞机 响应 问题 的 研究 不 至 于 过 分 
复杂 ,在 保证 模型 准确 度 的 基础 上 可 进行 适当 的 假设 。 一 般 来 说 ,大 气 率 流 是 一 个 随机 函 


< 密 的 关系 ,这 种 函数 关系 是 基于 大 量 的 况 
中 ,可 以 假设 :大 气 率 流 的 统计 特征 ( 即 平均 值 和 均 方差 ,以 及 相 


= 


量 和 统计 数据 来 进行 


KRANAR RO 既 不 随时 间 而 变 (认为 率 流 是 平稳 的 ) ,也 不 随 位 置 而 变 (认为 素 流 是 


均匀 的 ) ,在 这 种 假设 条 件 下 可 以 满足 对 飞机 品质 特性 分 析 的 需要 。 


[0125] 


[0126] 


数 型 纵向 相关 函数 为 :f (8) =e (基于 测量 


Dryden 率 流 模型 通过 有 限 差 分 方程 的 滤波 器 ,使 用 有 界 的 白 噪声 ,把 率 流 对 系统 
的 影响 施加 进 仿真 过 程 ,这 种 算法 首先 在 Military Specification MIL-F-8785C 和 
Military Handbook MIL-HDBK-1797 中 提出 。 


目前 常用 的 大 气 率 流 模型 是 德 莱 顿 基于 大 量 涡 


|= 


量 和 统计 数据 提出 的 ,其 对 应 的 指 


和 统计 数据 求 得 f O Mg (S) 后 ,然后 对 得 到 


的 结果 通过 Fourier 变 换 ,可 求 出 Dryden 模 型 的 纵向 和 横向 频谱 函数 如 下 : 


[0127] 


[0128] 


[0129] 


MORTA L 

z ]+(L,Q) 
= l (87) 

z 1+(L,Q) 
DQ) =0 l 

z 1+(L,02) 


lm ®(Q) = const 
Qo0 


im ®(O) x Q ° 


, 
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[0130] ”其 中 , 9 为 空间 频率 ,单位 为 弧度 / 米 ,L,,L,,L, 为 灌流 尺度 ;0,,0,,0, 为 灌流 强 度 。 
[0131] ”得 到 的 频谱 函数 是 基于 “平稳 ”和 “均匀 ”两 个 假设 条 件 得 到 的 ,这 就 会 使 得 所 得 到 
的 这 个 模型 在 无 穷 远 处 的 渐进 性 质 是 不 符合 实际 素 流 理论 的 ,但 这 并 不 影响 实际 的 工程 
应 用 ,该 Dryden 模 型 的 优点 在 于 :相对 来 说 频谱 的 形式 较为 简单 ,可 以 通过 常规 的 数学 方 
法 进行 处 理 ， 而 这 一 点 对 于 率 流 数值 仿真 是 非常 必要 的 

[0132] ”要 生成 大 气 率 流 数据 ,必须 要 使 用 成 型 滤波 器 进行 计算 ,以 伪 随 机 信号 作为 输入 ， 
而 且 成 型 滤波 器 的 设计 必须 依据 于 理论 频谱 函数 ,这 样 才能 够 保证 所 生成 的 大 气 亲 流 数 
据 频谱 特性 的 正确 性 .成 型 滤波 器 就 是 通过 滤波 环节 将 白 噪 声 转 化 为 期 望 的 有 色 噪 声 。 
[0133] ”根据 公式 : 


D, (0) -> HDI + 过 È arigos Go O] ag 
(i =1,2...,m) 
[0135] ”可 以 得 出 : 
[0136] 中 wo) 一 16GGjo) Ob (6) — (O 
[0137] ”由 此 可 见 , 只 要 把 给 定 的 输出 频谱 按照 上 式 分 解 ,就 可 以 得 到 成 型 滤波 器 的 传递 
函数 G (s) 。 低 空 扰 动 模型 使 用 的 是 Dryden 率 流 扰 动 模型 ,其 空间 域 频谱 表达 式 , 再 根据 方 
向 上 飞机 分 速度 V 和 频率 的 关系 : 
[0138] w=QV (41) 
[0139] F, o 为 时 间 频 率 , 得 到 模型 的 时 间 域 频谱 表达 式 为 : 
o (Q) 和 
z 1+(Lø@/V) 
2 L, 1+12(L,@/ V7 
" z [+4(L,o/Vy Í 
2 L, 1+12(L,@/ VF 
z +4, ovy} 
[0141] “对 上 述 各 式 进 行 分 解 , 则 可 得 到 为 了 产生 给 定 频 谱 @ ,(o ) 所 需要 的 成 型 滤波 器 
的 传递 函 数 G(s) .对 于 3 个 分 量 , 求 出 所 需 的 传递 函数 如 下 : 


[0134] 


V 
1 


l 
G S= t, 
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K 
G (s)= u 
ETA 
K, = Oy C3 u = 入 
\ ay V 
G,(s)= te 
(Tss +1) 
[0142] T š (43) 
K, = O, L aTa = h aTa L, 
\ aV V V 
G,(s)= Sd 
(Tas +1) 


K, = On L, /A = 243L, fp = 2L, 
V zV V V 


[0143] ”以 上 所 得 到 的 就 是 可 以 进行 大 气 率 流 数值 仿真 计算 所 需要 的 成 型 滤波 器 所 对 应 
的 传递 函 数 的 数学 表达 式 。 对 于 成 型 滤波 器 的 形式 还 不 能 直接 在 程序 中 运算 ,必须 要 由 滤 
波 器 所 对 应 的 传递 函数 的 形式 得 到 可 以 在 程序 中 以 数值 方法 进行 运算 的 递 推 公式 ,这 也 
就 是 滤波 器 在 程序 中 的 实现 方式 .对 于 低空 扰动 的 Dryden 率 流 扰动 模型 , 则 有 如 下 递 推 公 
I 


u: Ya = K,(xdt +T y,)/(T,, + dt) 
[0144] v: Yuz = K | rra), T dx +T +2T,dDy, Tay TE +Tadt+dr) A) 。 
wiy, = K [T dt di )x. -T dtx +(QOT2,+2T O dD)y. Thy, TZ +T dt + dt) 
[0145] (9) 离散 突 风 模 型 
[0146] “对 离散 突 风 来 讲 , 突 风 的 垂直 速度 变化 是 确定 的 ,属于 确定 性 突 风 , 也 就 是 说 其 突 
风速 度 不 随时 间 变 化 ,是 定常 的 。 通 常 我 们 用 到 的 离散 突 风 模型 有 两 种 :简单 点 的 锐 边 突 
风 模 型 以 及 复杂 点 的 1-cos 离 散 突 风 模型 .在 实际 应 用 中 ,由 于 1-cos 突 风 模 型 考虑 了 突 风 
场 的 空间 尺度 ， 在 评价 飞行 品质 和 飞 控 系 统 设计 时 ,广泛 采用 此 模型 
[0147] ”在 Military Specification MIL-F-8785C 中 提出 了 一 直 种 1-cos 形 状 的 阵风 模型 
(参见 图 8 所 示 ) 。 阵 风 可 单独 作用 于 每 一 个 轴 或 者 同时 作用 与 三 个 轴 。 在 模型 中 ,可 以 指定 
阵风 强 E (由 阵风 产生 的 风速 增加 ) ,阵风 长 度 (以 米 为 单位 的 阵风 长 度 ) ,以 及 阵风 开始 时 
间 。 离 散 阵 风 模 型 可 以 以 英尺 / 秒 和 米 / 秒 表示 风速 。 离 散 的 阵风 可 以 一 次 或 者 多 次 应 用 ， 
以 评估 飞机 多 大 风干 扰 的 响应 。 离 散 阵风 的 数学 表达 式 为 : 


0 (x<0) 
[0148] V. = wq i cos 75) (O<x<d,) (45) 
y, (> 不 ] 


[9149] 其 中 ,V 为 阵风 强度 ,x 为 阵风 长 度 ,V ,为 作用 在 体 轴 轴 线 方向 上 的 最 终 风速 。 
[0150] (D 大 气 参数 模型 

[0151] ”1976 年 美国 Committee on Extension to the Standard Atmosphere (COESA) 机 
构 发 布 了 美国 标准 大 气 参 数 , 这 一 参数 仅 与 飞行 器 所 在 高 度 有 关 , 包 含 了 大 气 绝 对 温度 、 
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大 气压 力 、 大 气 密 度 和 声速 ,在 32000 米 以 下 ,美国 标准 大 气 与 国际 民航 组 织 定 义 的 标准 大 
气相 同 。 参 见 图 9 所 示 。 

[0152] ”飞行 高 度 大 于 84852 米 时 ,温度 值 是 由 线性 外 推 得 到 ,压力 值 是 由 对 数 外 推 得 到 ， 
密度 和 声速 的 值 则 由 理想 气体 方程 得 到 。 

[0153] (5) 重力 加 速度 模型 

[0154] WGS-84 坐 标 系 (World Geodetic System 一 1984Coordinate System) 是 一 种 国际 
上 采 用 的 地 心 坐标 系 。 原 点 是 地 球 的 质心 ,空间 直角 坐标 系 的 Z 轴 指向 BIH(1984.0) 定 义 的 
地 极 (CTP) 方 向 , 即 国际 协议 原点 CI0, 它 由 IAU 和 IUGG 共 同 推 荐 .Xx 轴 指向 BIH 定 义 的 零度 子 
午 面 和 CTP 赤 道 的 交点 ,Y 轴 和 Z,X 轴 构成 右手 坐标 系 。.WGS-84 顶 球 采 用 国际 大 地 测量 与 地 球 
物理 联合 会 第 17 届 大 会 测量 常数 推荐 值 ,采用 的 两 个 常用 基本 几何 参数 。 如 图 10 所 示 。 
[0155] ”WGS-84 是 修正 NSWC9Z-2 参 考 系 的 原点 和 尺度 变化 ,并 旋转 其 参考 子午 面 与 BIH 定 
义 的 零 度 子午 面 一 致 而 得 到 的 一 个 新 参考 系 ,WGS--84 坐 标 系 的 原点 在 地 球 质心 ,Z 轴 指向 
BIH1984.0 定义 的 协定 地 球 极 (CTP) 方 向 ,Xx 轴 指 向 BIH1984.0 的 零度 子午 面 和 CTP 赤 道 的 交 
点 ,Y 轴 和 Z、X 轴 构成 右手 坐标 系 。 它 是 一 个 地 固 (地 心 固 连 ) 坐标 系 , 有 公式 : 

[0156] ”长 半径 :a 二 6378137 土 2 (m) ; 

[0157] “地球 引力 常数 和 地 球 质 量 的 乘积 :GM 王 3986005X108ms-2 士 0.6X108mss *; 

[0158] ”正常 化 二 阶 带 谐 系 数 :C,, 二 -484.16685X10-6 土 1.3X10-9; 

[0159] ”地 球 重力 场 二 阶 带 球 谐 系 数 :J, 二 108263 X10-8 

[0160] ”地 球 自转 角速度 : o 二 7292115X10-1llrads ' 填 0.150X10-1llrads ' 

[0161] ” 扁 率 f= 二 0.003352810664 

[0162] ”在 World Geodetic System(WGS84) 中 ,根据 经 度 、 纬 度 和 高 度 ,在 地 球 的 等 位 椭 球 
上 确定 重力 加 速度 。 

[0163] ”重力 加 速度 和 飞机 所 在 高 度 的 关系 ,参见 图 11 所 示 。 

[0164] (三 ) 对 飞机 所 受气 动力 以 及 飞机 气动 力 的 弹性 影响 建 模 

[0165] ”准确 的 气动 力 分布 是 进行 准确 静 气 动弹 性 特性 分 析 和 飞机 型 架设 计 的 基础 。 飞 机 
不 同 设计 阶段 有 不 同 的 精度 要 求 , 需 要 采用 不 同 的 气动 分 析 方 法 ,在 飞机 设计 初期 ,快速 
高 效 地 形成 设 计 方 案 , 确 定 气动 主要 参数 ;飞机 详细 设计 阶段 需要 计算 准确 的 弹性 气动 
力 ,为 相关 专业 提供 气动 数据 输入 。 因 此 ,采用 满足 飞机 所 在 设计 阶段 精度 需求 的 弹性 气 
动 流 场 求解 方法 十 分 必要 。 
[0166] ”气动 系数 计算 根据 飞行 器 的 基本 构 型 ,分 别 计算 气流 坐标 系 下 纵向 基本 气动 系 
数 、 横 航向 基本 气动 系数 .升降 能 产生 的 气动 系数 、 副 翼 产 生 的 气动 系数 方向 能 产生 的 气 
JRA RE 产生 的 气动 系数 、 扰 流 板 产生 的 气动 系数 、 动 导数 产生 的 气动 系数 ,然后 计算 
各 部 分 的 气动 系 数 之 和 作为 飞机 整体 的 气动 系数 。 

[0167] ”计算 出 气流 坐标 系 下 的 气动 系数 后 ,然后 根据 气流 角 ( 迎 角 、 侧 滑 角 ) 将 气动 系数 
转化 到 机 体 坐 标 系 下 。 

[0168] ”根据 飞行 器 当前 的 动 压 ( 跟 飞行 速度 、 高 度 相关 ) 、 重 心 位 置 计算 气动 力 和 力矩 沿 
机 体 轴 的 分 量 。 
[0169] ”气动 力 系 数 数据 库 

[0170] ”在 本 案例 中 , 先 建立 相关 参数 关于 气动 系数 及 气动 导数 的 多 维 数据 库 , 利 用 参数 
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辨识 校 核 风 洞 实验 数据 ,相关 结果 与 实验 值 的 对 比如 下 ,其 中 ,升力 计算 结果 在 实验 点 满 
足 5% 的 误差 要 求 。 阻 力 系数 基本 满足 误差 范围 要 求 ; 除 了 在 小 攻 角 (0 "附近 ) 由 于 绝对 值 
很 小 导致 相对 误差 略 有 超出 。 升 力 与 迎 角 的 关系 曲线 如 图 12 所 示 。 阻 力 与 迎 角 的 关系 曲线 
如 图 13 所 示 。 
[0171] ”在 仿真 运行 时 ,可 以 对 所 有 气动 参数 进行 插值 ,进而 直接 得 到 得 到 当前 迎 角 、 操 纵 
面 下 机 如 的 气动 参数 ,从 而 输入 下 一 模块 进行 力 和 力矩 的 设计 ,最 后 与 六 上 自由 度 非 线 性 运 
动 模块 相 耦 合 , 构 建 信 息 反 馈 。 
[0172] ”地 面 效 应 对 气动 力 的 影响 
[0173] ”从 简化 的 空气 动力 学 理论 看 ,地 面 对 于 单独 机 经 或 尾 咒 相当 于 一 个 镜 象 涡 系 ) 对 
此 简化 模 型 的 理论 分 析 和 相应 的 实验 研究 早 在 40-50 年 代 即 已 完成 ,然而 作为 机 经 、 机 身 
和 平 尾 的 组 合 体 , 并 考虑 到 地 面 对 大 偏 度 襟 细 和 干扰 流 场 及 配 平 力矩 的 综合 影响 ,情况 就 
杂 得 多 ,目前 尚 无 完整 和 可 靠 的 理论 计算 方法 ,由 于 雷诺 数 及 其 它 实 验 条 件 (如 活动 地 
板 和 配 平 舵 偏 角 ) 的 限制 ， 风 洞 实验 的 精度 也 受到 限制 。 
[0174] ”这 里 由 镜 象 涡 理论 配合 经 验 修正 的 综合 方法 出 发 ,考虑 到 飞行 力学 中 的 力矩 配 平 
关系 , 形 成 了 计 及 地 效 时 全 机 的 附加 升力 、 阻 力 、 迎 角 和 下 洗 角 的 工程 估算 方法 ,在 给 定 升 
力 系 数 CL 和 离 地 高 度 h (平均 气动 弦 1/4 处 距 地 高 度 ) 后 ,由 尾 涡 引起 的 迎 角 减少 量 (弧度 ) 
和 附着 涡 引 起 的 升力 系数 变化 量 各 为 : 


GC 
[0175] ôq = -一 于 
TÅ (46) 
1 
[0176] óC, = 一 一 一 | (47) 
x 0 rey | 
[0177] ”其 中 ,各 参量 的 计算 公式 为 : 
0.788 
[0178] c = exp -2.43 25 | (48) 
[0179] y= [2] -2 (49) 
b b 
he 
[0180] C50) 


g 8z| (7/2) +(1/64) | 


[0181] JErh,h 234 05J9Z1/41byB Hb) J b WRR AFADI N 3š G HASE 
1/452 A BJ JIE #280Cm1 / 407E EPE A 10 t se GERTA JJ kE # 280 m1/4; 


[0182] N =1l+(N'-1) I+ (G, -cu (51) 
L 


[0183] N? AERE FE | eW samj , UEAN : 


2 8 
108 48 | LI (52) 
h h p 


[0184]  N'=1+0.00239 
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[0185] 地面 效 应 影响 下 的 升力 曲线 如 图 14 所 示 。 

[0186] ”考察 地 面 效 应 导致 的 飞机 升力 系数 增 量 ,与 高 度 的 关系 如 下 图 15 所 示 。 

[0187] ”气动 弹性 模块 设计 

[0188] ”为 了 减少 计算 量 ,考虑 巡航 飞行 状态 下 的 飞行 力学 配 平 ,在 气动 弹性 的 计算 中 , 含 
去 了 对 计算 结果 影响 不 大 的 部 分 ,首先 根据 当前 的 气动 参数 ,计算 机 可 上 的 气动 载荷 ,再 
将 集中 载荷 根据 经 验 公式 分 布 在 机 辟 表 面 , 从 而 计算 机 翼 的 二 维 气动 -结构 耦合 的 结果 ， 
再 将 机 经 变形 量 输 回 到 气动 力 模块 ,得 到 当前 状态 下 的 弹性 影响 的 气动 力 。 

[0189] “大 展 弦 比 的 升力 面 EIE, 后 掠 辟 ) , 当 其 结构 的 弹性 特性 可 用 工程 梁 理论 表达 
时 , 则 是 一 维 静 气动 弹性 问题 .其 定常 气动 力 可 以 用 修正 片 条 理论 ,如 考虑 沿 翼 展 下 洗 流 
变化 , 则 可 以 用 L.Prandtl1 和 J.Weissinger 的 升力 线 理论 ,当然 把 工程 问题 简化 为 一 维 问 
题 处 理 , 无 论 气动 力 和 弹性 力 (变形 ) , 均 做 一 定 简化 假设 ,为 求解 更 精确 ,大 展 弦 比 升力 面 
也 可 以 作为 二 维 问题 处 理 。 
[0190] ” 解 一 维 静 气动 弹性 问题 可 用 微分 方程 和 积分 方程 ,由 于 微分 方程 时 其 气动 力 只 能 
用 片 条 理 论 。 而 积分 方程 可 以 用 片 条 或 升力 线 理论 , 故 工 程 上 认为 用 积分 方程 矩阵 数值 解 
法 更 加 合理 。 
[0191] ”根据 图 16 所 示 的 机 辟 几 何 模型 ,经 推导 ,积分 方程 为 : 
[0192] ”对 直 机 经: 


piny 


Í s e S 
0(y)= af” ined c; Te n+ f0) 
0 


[0193] (53) 


f OC ¿ 
f(y)= afela] Cuac +e yan |C (amy medn 
0 0 


[0194] X EHR, ARHAN : 


0(y)= dc %(y,n)+ :dc T sh n+ f(y) 
K e oð 


[0195] ; m ; i (54) 
广 O)=d| co， as + “espa = fle On), n) nmgan 
0 0 
0 
a(y)=a,(y)+0(y)+ 6 Ü) 
[0196] I: ac (55) 
C,(y)=C;(y)+C;(y)+ó 全 (y) 


[o197] 后 掠 轩 与 直 机 村 不 同 ,弯曲 变形 可 以 引起 顺 流向 气流 迎 角 的 变化 ,所 以 C”(y ,由 和 
C G 略 均 分 别 与 全 式 梁 的 弯曲 、 扭 转 刚度 有 关上 面 各 式 中 

[0198] 9: 远方 来 流速 压 g= 了 PP2;6 :操纵 面 偏转 角 ,rad:0 0) :y 处 剖面 顺 航向 弹性 变形 
farad; C" (yn) :结构 影响 系数 ,剖面 作用 单位 力矩 0 方向 ) 引起 y 剖 面 的 弹性 变形 角 ( 
方向 ) ; C Gy,m :结构 影响 系数 ,n 剖 面 作用 单位 力 (z 方 向 ) 引起 y 剖 面 的 弹性 变形 角 (0 广 
向 ) ，a cvevdimvne sg: 分 别 指 a 站 位 的 几何 迎 角 、 纺 长 ` 气 动 中 心 至 刚 轴 的 距离 ,可 段 重 必 
至 刚 轴 的 距离 ,单位 展 长 (y 方 向 ) 质量 ,2 向 过 载重 力 加 速度 : CI, CEC, ic: 气动 力 系数 ,下 
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标 L 表 示 升 力 ,M 表 示 力 矩 ,AC 表 示 对 气动 中 心 , 上 标 r 表 示 由 刚性 、e 表 示 由 弹性 变形 引起 
的 。 
[0199] ”如 定义 新 的 柔 度 影响 系数 : 

C'(p,n)= C2(y,n)+— C%(y,n) 
[0200] 了 (56) 

cm) C” (y,n)+ 7 C“ (n) 


[0201] D ELR JA BR KOKERE H i T FI Bi — 088 Jr Bak 2: 


{}= dllcc +E + 一 中 +{f} 

[0202] ge (57) 
Vyd hcnt al- ane) 

[E] = |c; lo, fe] 


[9203 stt: [F] leo, fe] 
[G] = [c; Ë. ja! (58) 
[0204] ”其 中 区 | 为 数值 积分 加 权 和 矩阵, 如 为 等 间距 间隔 ,可 用 Simpson 或 梯形 法 则 ,然而 用 


升力 线 理 论 时 , 展 向 均 非 等 间隔 , 求 数值 积分 有 较 高 难度 ,可 用 Multhopp 求 积 方程 , 按 图 17 
把 区 间 (a,b) 划分 出 Multhopp 站 位 , 则 按 Multhopp 求 积 方程 : 


f _ z(b-a) x . ! 
[201 froe ir 227/)sn $, (59) 
[0206] PA Jz2Multhopp12 k P šë E HIRI 


ZEC 


[0208] 即 可 计算 得 到 |w ,| 矩阵 。 


[0209] ”根据 以 上 公式 , 迎 角 2 ,马赫 数 为 0.8 时 飞机 升力 系数 治 展 向 分 布 如 图 18 所 示 。 可 
以 看 出 , 当 飞 机 机 可 为 弹性 体 时 ,机 可 治 展 问 发 生 扭 转 , 飞 机 升力 系数 较 刚 性 状态 下 略 有 
偏 大 。 

[0210] (四) 将 上 述 三 个 模型 连 成 通路 ,建立 模拟 系统 ,最 终 得 到 弹性 飞机 的 阵风 相应 。 
[0211] ”基于 所 建立 的 飞机 本 体 、 环 境 和 气动 力 模 型 ,将 以 上 模型 按 模 拟 过 程 中 的 物理 规 
则 连接 起 来 ,首先 ,在 每 次 模拟 中 ,根据 六 自由 度 模块 得 到 的 飞机 飞行 状态 量 , 经 环境 模块 
计算 得 到 飞 机 当地 所 处 的 风速 ,温度 、 重 力 加 速度 等 物理 变量 ;将 这 些 物理 变量 与 飞机 飞 
行 状态 量 结 合 ， 得 到 气流 坐标 系 中 的 飞机 速 压 、 迎 角 、 侧 滑 角 等 信息 ;把 这 些 信息 放 入 气动 
力 数据 库 进 行 插值 ， 得 到 飞机 的 气动 系数 ;在 地 面 效应 和 气动 弹性 模块 的 作用 下 ,气动 系 
数 的 值得 到 修正 ;最 后 和 重力 一 起 ,计算 飞机 飞行 的 合力 以 及 合力 矩 , 传 回 六 自由 度 模 块 ， 


(60) 


f 
[0207] | f (yy = 
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计算 下 一 时 刻 的 飞机 飞行 状 态 ,形成 通路 ,如 图 19 所 示 为 Simulink 模 块 的 连接 关系 ， 

本 发 明 的 特点 在 于 把 飞行 动力 学 、 飞 机 弹性 变形 和 空气 动力 三 者 结合 ,分 别 建立 
模型 , 满 足 实 时 仿真 要 求 ,建立 统一 的 弹性 飞机 飞行 动力 学 仿真 平台 ,针对 任意 弹性 飞机 
和 任意 飞行 状 态 ,得 到 对 应 的 全 波长 和 半 波 长 阵风 相应 ,为 飞机 的 柔性 设计 提供 了 统一 的 
仿真 平台 ,减少 设 计 投 入 ,缩短 设计 周期 。 


[0212] 


附 图 说 明 

[0213] ”图 1 为 气动 伺服 弹性 力学 三 角形 图 示 。 

[0214] ”图 ?为 Collar 气 动弹 性 三 角形 图 示 。 

[0215] ”图 3 为 飞机 运动 坐标 系 图 示 。 

[0216] ”图 4 为 风 模 块 和 率 流 模块 。 

[0217] ”图 5 为 风 矢 量 在 垂直 方向 上 的 变化 。 

[0218] ”图 6 为 切 变 风 在 竖 直 方向 上 的 分 布 。 

[0219] ”图 7 为 率 流 强度 的 均 方 差 和 高 度 的 关系 。 

[0220] ”图 8 为 1-cos 突 风 模 型 。 

[0221] ”图 9 为 美国 标准 大 气 1976。 

[0222] ”图 10 为 WGS-84 坐 标 系 图 示 。 

[0223] ”图 11 为 重力 加 速度 和 飞机 所 在 高 度 的 关系 。 

[0224] ”图 12 为 升力 与 迎 角 的 关系 曲线 。 

[0225] ”图 13 为 阻力 与 迎 角 的 关系 曲线 。 

[0226] ”图 14 为 地 面 效应 影响 下 的 升力 曲线 。 

[0227] ”图 15 为 地 面 效 应 导致 升力 系数 增 量 与 高 度 关系 。 

[0228] KIOJ JUTRA. 

[0229] ”图 17 为 把 区 间 (a,b) 划分 出 Multhopp 站 位 。 

[0230] ”图 18 为 弹性 飞机 与 刚性 飞机 升力 系数 沿 展 向 分 布 对 比 。 

[0231] ”图 19 为 搭建 的 整体 仿真 系统 

[0232] ”图 20 为 某 小 型 通用 飞机 几何 模型 。 

[0233] ”图 21 为 某 小 型 通用 飞机 中 作用 在 副 辟 重心 上 的 载荷 随时 间 变 化 曲线 。 
[0234] ”图 22 为 某 小 型 通用 飞机 中 作用 在 飞机 重心 上 的 载荷 随时 间 变 化 曲线 。 
[0235] ”图 23 为 某 大 型 客机 几何 模型 。 

[0236] ”图 24 为 某 大 型 客机 中 辟 梢 位 移 随时 间 变 化 曲线 。 

[0237] ”图 25 为 某 大 型 客机 迎 角 随 时 间 变 化 曲线 。 

具体 实施 方式 

[0238] 下面 对 该 阵风 相应 分 析 方 法 进行 案例 验证 。 

[0239] (一 ) 某 小 型 通用 飞机 弹性 阵风 相应 的 模拟 验证 

[0240] “在 某 小 型 通用 飞机 的 弹性 阵风 相应 的 分 析 中 ,飞机 几何 外 形 如 图 20 所 示 , 机 可 为 


平 直 下 单 既 ,由 于 飞机 较 小 ,机 如 、 尾 吧 的 弹性 变形 相差 不 大 ,因此 需要 同时 对 机 器 、 尾 贵 
进行 弹性 建 模 ,考虑 弹性 影响 ,在 全 波长 离散 阵风 的 激励 下 ,飞机 状态 经 过 1 秒 左右 的 震荡 
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逐渐 回 到 稳 定位 置 ,其 中 ,作用 在 副 妈 重心 上 的 载荷 随时 间 变 化 曲线 以 及 作用 在 飞机 重心 
上 的 载荷 随时 间 变化 曲线 如 图 21、 图 22 所 示 。 
[0241] (二 ) 某 大 型 客机 弹性 阵风 相应 模拟 验证 
[0242] ”在 某 大 型 客机 的 弹性 阵风 相应 的 分 析 中 ,飞机 几何 外 形 如 图 23 所 示 , 机 骂 为 后 拯 
下 单 细 ， 由 于 飞机 较 大 ,机 经 的 弹性 变形 显著 大 于 全 机 其 他 位 置 , 因 此 仪 需 要 对 机 如 进 行 
弹性 建 模 , 考 虑 其 弹性 影响 ,在 半 波 长 离散 阵风 的 激励 下 ,飞机 状态 经 过 8 秒 左右 的 震荡 未 
渐 回 到 稳定 位 置 ， 其 中 ,经 梢 位 移 随时 间 变 化 曲线 以 及 飞机 迎 角 随时 间 变 化 曲线 如 图 24、 
图 25 所 示 。 
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